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ABSTRACT 


A  two-dimensional  ray  tracing  analysis  for  the  calculation  of 
radome  horesight  error  and  antenna  pattern  distortion  is  presented  here. 
£}mphasis  has  been  placed  on  the  development  of  a  method  having  con¬ 
siderable  flexibility,  so  as  to  enable  application  of  the  method  to  a  wide 
range  of  antenna- radome  problems,  and  on  relative  ease  of  calculation, 
so  as  to  niimmize  calculation  time.  Several  example  problems  are  cal¬ 
culated  to  demonstrate  the  usefulness  of  the  approach.  Comparisons 
between  calculations  and  measurements  have  been  included  whenever 
measured  dataware  available.  Instructions  for  use  of  this  completely 
computerized  method  are  included  along  with  several  tables  describing 
variables  and  the  complete  computer  program  with  necessary  sub¬ 
routines.  Programs  are  written  in  Fortran  IV  language  suitable  for 
use  on  the  OSU  version  of  the  IBM  system  360/ 75  (some  minor  changes 
may  be  required  for  use  on  other  360/75  installations). 
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A  TWO-DIMENSIONAL  RAY- TRACING  METHOD  FOR  THE 
CALCULATION  OF  RADOME  BORESIGHT  ERROR 
AND  ANTENNA  PATTERN  DISTORTION 


I.  INTRODUCTION 

Streamlined  radomes  for  aircraft  and  missile  guidance  systems 
must  be  carefully  designed  for  high  transmission  efficiency  and  minimum 
boresight  error.  Since  the  usual  antenna- radome  system  is  large  in 
terms  of  wavelengths,  exact  methods  foi-  the  calculation  of  radome  errors, 
such  as  the  integral  equation  methods  of  Van  Doeren^  and  Hahn,*  prove 
to  be  difficult  to  apply.  Frequently  these  methods  can  only  be  applied 
to  a  small  portion  of  the  radome  such  as  the  vertex  region.  Therefore 
approximate  methods  continue  to  be  useful  in  radome  analysis. 

This  report  presents  a  two-dimensional  approximate  method  for 
calculating  radome  boresight  error  ar.d  antenna  pattern  distortion.  A 
ray  analysis  is  used  to  determine  the  effects  of  the  radome  on  the  antenna. 
These  effects  are  used  to  modify  the  source  aperture  distribution  which 
is  nimierically  integrated  to  determine  the  far- zone  field  pattern  of  the 
antenna- radome  system.  The  Ohio  State  University- IBM  System  360/75 
high  speed  digital  computer  is  used  for  all  calculations.  The  calculated 
results  agree  reasonably  well  with  experimental  data  and  require  little 
computer  time.  Several  calculations  of  typical  radome  design  problems 
are  discussed. 


II.  THE  BASIC  METHOD 

A.  General  Considerations 


The  analysis  is  based  upon  a  two-dimensional  model  of  the  antenna- 
radome  system  as  shown  in  Fig.  1.  The  radome  is  represented  by  its 
cross-section  and  the  source  antenna  is  represented  by  a  one-dimensional 
aperture  having  a  known  amplitude  and  phase  distribution.  Rays  are 
traced  from  the  aperture  to  the  radome  wall  to  determine  angles  of 
incidence  to  be  used  in  calculating  radome  effects.  The  radome  is 
approximated  by  a  plane  multilayer  oriented  at  the  calculated  angle  of 
incidence  at  each  ray  intersection.  The  plane  wave,  plane-sheet  trans¬ 
mission  coefficient  and  inseition  phase  delay  are  calculated  for  each 
ray.  These  values  are  used  to  modify  the  original  source  distribution 
function  such  that  a  reconstructed  aperture  distribution  is  obtained  which 
includes  the  radome  effects.  This  distribution  function  is  then  numeri¬ 
cally  integrated  by  high-speed  digital  computer  to  determine  the  approxi¬ 
mate  far-field  pattern  of  the  antenna- radome  system  which  is  compared 
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Fig.  1.  Two-dimensional  model  for  antenna- radome  system. 


to  the  pattern  obtained  without  the  radome  to  determine  pattern  dis¬ 
tortion  and  boresight  error. 

®  •  Ray  Tracing 

In  the  usual  antenna- radome  system  the  antenna  aperture  plane 
is  displaced  by  some  distance  da  from  the  gimbaling  axes  of  the  antenna. 
When  the  antenna  is  scanned,  description  of  the  aperture  plane  becomes 
difficult  in  a  fixed  coordinate  system.  For  this  reason  two  coordinate 
systems  are  used  to  describe  the  antenna- radome  geometry,  as  shown 
in  Fig.  2,  The  radome  is  described  in  a  fixed  (x,  y)  frame  which  has  its 
axes  centered  on  the  antenna  gimbal  axis.  The  antenna  aperture  is 
described  in  {x',y')  frame  which  rotates  about  the  antenna  gimbal  axis 
with  the  angle  of  rotation  corresponding  to  some  look  angle  <J)L,  Points 
in  the  (x',y')  system  are  related  '  points  in  the  (x,y)  system  by  the 
following  transformation: 

(1)  ^x\_/cos4>l  -  sin<))L 

Vy;  \sin<j>L  cosfpL,;  ^y'j 


The  radome  is  asstimed  to  be  constructed  of  n  geometry  sections 
which  can  be  described  by  the  following  general  secondrordei  equation: 


(2)  F{x,  y)  =  anx^  +  bny^  +  Cnxy  +  dnX  +  enY  +  fn  =  0 


where  (x,  y)  are  the  coordinates  of  Fig.  2  and  a^.  •  .fn  are  a  set  of  geo¬ 
metrical  constants  which  define  the  n-th  radome  section.  A  set  of  m 
equally- spaced  rays  from  the  antenna  aperture  to  the  radome  inner  wall 
ave  selected  to  represent  the  problem.  A  ray  drawn  from  a  point  (xa,  ya) 
on  the  aperture  plane  to  the  radome  wall  is  described  by  the  point-slope 
form  as: 


(3)  y  -  ya  =  mR(x  -  xa) 


where  mR  is  the  slope  of  a  ray  in  the  (x,  y)  frame.  Th  antenna  points 
to  be  used  vare  determined  in  the  (x',y')  frame  by 


(4) 
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Fig.  2.  Coordinate  system  used  to  define  the 
antenna- radome  geometry. 


(5)  y^  =  — (2m  -  2  -  NR) 

2  NR 


where 


A  is  the  total  aperture  length 
NR  is  the  number  of  rays  to  be  used 
m  is  the  index  of  a  particular  ray 
da  is  the  perpendicular  distance  from  the  origin  to 
the  aperture  plane  as  in  Fig.  2. 
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The  set  of  m  points  determined  by  Eqs .  (4)  and  (5)  are  transformed  by 
Eq.  (!)  to  the  set  of  {xai  ya)  points  to  be  used  in  Eq.  (3).  Substituting 
Eq.  (3)  into  Eq.  (2)  we  get  the  following  quadratic  in  x: 


(6) 


X 


m 


(-2bmR  ^a  2bmj^ya  -  cmRXa  +  cya  +  d  +  emR) 
(a  +  bmR'^  +  cmR) 


(-2bmRyaXa  +  bmR  Xa^  +by^  -  emRX^  +  eya  +  f) 

+ - ^ -  =  0 

(a  +  bmj^^  +  croR) 

The  solution  of  Eq.  (6)  gives  the  x-coordinate  of  the  point  of  intersection 
of  the  m-th  ray  and  the  radome.  Since  Eq.  (6)  is  of  the  form: 


(7)  +2  Bxjn  +  C  =  0 


the  solution  of  Fn.  (6)  can  be  written  as 

(8)  x„i  =  -  -  C 

where 


2B  is  the  coefficient  of  the  linear  term  in  Eq.  (6) 

C  is  the  constant  term  in  Eq.  (6)  . 

From  the  geometry  of  the  system  it  is  seen  that  the  positive  square 
root  is  selected  in  Eq.  (8)  to  give  the  proper  point  of  intersection.  This 
value  of  Xm  is  substituted  into  Eq.  (3)  to  obt  .n  the  y-coordinate  of  the 
intersection  point: 


(9)  Yn,  =  mR(xn,  -  Xa)  +  Ya  • 


The  derivative  of  Eq.  (2)  evaluated  at  {xin>ym)  gives  the  slope  of 
the  tangent  to  the  radome  surface  at  the  m-th  intersection  point; 


(10) 


(2  a  Xj;yj  +  C  Yi^yj  +  d) 
mT  = - - - - - 
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Provided  that  neither  the  tangent  line  to  the  radome  nor  the  ray  is 
parallel  to  the  y-axis  and  that  the  two  lines  are  not  perpendicular,  the 
angle  of  intersection  of  the  two  lines  is: 


(H) 


eP  -  Tan-» 


-  nrirp 

1  +mRm'p 


which  is  the  complement  of  the  angle  of  incidence  of  the  m-th  ray  and 
the  radome  inner  wall. 


(12)  eP  =  Z  -  em 

'  I  2 


is  the  angle  of  incidence.  As  will  be  discussed  later  the  average  angle 
of  incidence  for  two  adjacent  rays  will  be  used  in  farther  calculations: 


(13)  6^=  (ef  +  0f^^)/2 


The  two  exceptions  to  Eq.  (H)  mentioned  above  are  treated  specifi¬ 
cally  in  the  computer  program.  The  constants  a  through  f  in  Eq.  (1) 
depend  upon  the  specific  geometry  of  the  radome.  Logic  statements  in 
the  program  assure  that  the  ray  intersection  is  calculated  in  the  proper 
geometrical  section.  The  angle  of  incidence  calculated  in  Eq.  (13)  is 
stored  in  an  m  X  n  array  indicating  that  the  m-th  ray  ic  used  with  the 
n-th  set  of  geometrical  and  electrical  constants.  The  an...fn  constants 
and  the  associated  n  geometry  boundaries  are  usually  calculated  in  the 
program,  however,  for  specialized  cases  they  may  be  read  in  directly. 


C.  Pattern  Calculation 

The  basic  calculation  is  tha*  of  a  se-tion  of  the  far  zone  field 
pattern  of  the  antenna  radiating  in  the  presence  of  the  radome.  The 
angular  range  over  which  the  pattern  .s  calculated  varies  from  one 
degree  about  the  antenna  look  angle  for  boresight  error  calculation  and 
from  10  to  90  degrees  about  the  look  angle  for  pattern  distortion  calcula¬ 
tion.  Within  the  one  degr-  e  interval  used  for  boresight  calculation  only 
a  few  discrete  pointc  are  calculated. 

Ihe  far-zene  field  pattern  for  the  one-dimensional  source  repr»-- 
sentation  shown  in  Fig.  1  is  given  by; 
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(14) 


E(4.)=  \  F(y)  ei^l' dy 


where: 

F(y)  is  the  amplitude  distribution  function 
‘l>(y)  is  the  phase  distribution  function 
4)  is  the  pattern  angle 
L  is  the  length  of  the  aperture . 

In  general  F(y)  and  4>{y)  are  arbitrary  func’ions  such  that  the 
evaluation  of  the  integral  requires  numerical  methods .  These  functions 
are  determined  by  the  given  source  distribution  functions  and  modified 
later  to  account  for  the  presence  of  the  radome.  Rays  are  traced  from 
the  aperture  plane  to  the  radome  inner  wall  where  they  are  modified  by 
the  plane  wave,  pla!ne-sheet  transmission  coefficient  (|t!“)  and  insertion 
phase  delay  (IPD),  to  a  new  aperture  plane  immediately  outside  the  radome. 
Here  a  "reconstructed”  aperture  is  defined  which  determines  the  far-field 
of  the  antenna- radome  system  according  to  Eq.  (H).  A  few  comments 
on  ray  tracing  follow. 

The  usual  ray  analysis  uses  a  set  of  n  equally- spaced  rays.  As 
this  n  is  increased  the  predicted  result  varies  up  to  some  value  of  n 
where  further  addition  of  rays  no  longer  changes  the  answer.  This 
answer  is  not  necessarily  the  correct  answer  but  merely  the  best  answer 
that  the  ray  tracing  solution  can  predict.  This  n  required  for  a  con¬ 
vergent  answer  using  equally  spaced  rays  frequently  becomes  quite  large, 
typically  500  rays  for  a  lOX  aperture  ana  a  streamlined  radome.  Evalu¬ 
ation  of  Eq.  (14)  using  a  large  n  consumes  =»n  excessive  amount  of  com¬ 
puter  time  which  is  undesirable.  An  alternative  to  this  approach  is  to 
use  a  set  of  fewer  unevenly- spaced  weighted  rays^  to  analyze  the  problem 
such  as  the  set  shown  in  Fig.  1.  Since  the  radome  effects,  |tI^  and  IPD, 
are  strongly  dependent  on  incidence  angle,  close  spacing  of  the  rays  is 
required  only  if  the  radome  curvature  is  changing  rapidly.  More  widely 
spaced  rays  can  be  used  in  regions  where  the  curvature  variation  is 
slight.  The  most  efficient  ray  analysis  uses  the  fewest  nvunber  of  rays 
required  to  obtain  the  convergent  answer.  The  numerical  treatment 
which  will  be  applied  tc  the  reconstructed  aperture  is  equivalent  to  per¬ 
forming  such  a  weighted  ray  tracing.  This  treatment  is  described  below. 

The  source  antenna  has  associated  with  it  a  large  number  of  equally- 
spaced  rays,  say  500.  A  subaperture  cf  the  source  is  defined  as  that 
section  of  the  source  aperture  between  two  successive  rays  and  has 
associated  v^ith  it  a  ray  emanating  from  its  center  which  intersects  the 
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radome  inner  wall  at  angle  6^  of  Eq.  (13).  The  local  amplitr.de  and 
phase  associated  with  this  subaperture  is  calculated  from  the  known 
source  distribution  function.  In  the  case  of  a  circular  aperture,  as 
shown  in  Fig.  3,  the  equivalent  one  dimensional  source  must  be  tapered 
by  the  factor 


(15) 


Ao(y)  =  cos 


This  factor  takes  into  account  the  effective  power  radiated  from  each 
segment  represented  by  the  chord  length  at  the  coordinate  y,  as  in 
Fig.  3.  If  the  circular  aperture  itself  has  an  amplitude  taper,  Fir),  the 
equivalent  one -dimensional  aperture  taper  required  is: 


(16)  A{y)  =  F(y)  Ac{y) 

where  Ao(y)  is  found  from  Eq.  (15). 


Fig.  3.  Amplitude  la^^cr  of  a  one-dimensional  aperture  for 
equivalence  to  a  two-dimensional  circular  aperture. 
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This  subaperture,  according  to  con/entional  ray  optics,  illumi¬ 
nates  only  the  small  subsection  of  the  raUome  wall  lying  between  its 
fwo  defining  rays  which  is  approximated  as  a  plane  sheet  oriented  at 
©A*  Plane  wave,  plane-sheet  |t1^  and  IPD  are  calculated  using  the 
method  of  Richmond*  for  the  rays  associated  with  each  subaperture. 

The  local  subaperture  field  distribution  is  modified  by  the  local  [t}^ 
and  IPD.  The  reconstructed  aperture  is  thus  completed  specifying  the 
integrand  function.  F(y)  and  <j>{y)  of  Eq.  (14). 

The  numerical  treatment  of  the  aperture  integral  involves  bresiking 
the  integral  down  into  several  sections,  or  subaperturas,  as  determined 
by  the  rate  of  change  of  the  integrand,  integrating  over  these  subaper¬ 
tures,  and  summing  the  integrals.  Equation  (17)  specifies  the  calcu¬ 
lation. 

N 

Fn(y)  dy 


where: 

l^n  •  length  of  the  n-th  subaperture 
4'n  =  phase  of  the  n-t?i  subaperture 
Fn  =  amplitude  of  the  n-th  subaperture 
<{)  =  pattern  angle  . 

The  process  of  determining  the  apertu'»'e  subdivision  is  as  follows. 
Fixed  amplitude  and  phase  deviations  are  specified,  usually  0.05  to 
0,10  and  2  to  3  degrees  respectively.  The  length,  amplitude,  and  phase 
of  the  first  subaperture  are  determined  by  scanning  from  the  center  of 
the  reconstructed  aperture,  point  by  point  towards  the  positive  endpoint 
of  the  aperture,  until  either  of  the  fixed  deviations  occurs.  At  this 
poiiit  the  first  subaperture  boundary  is  defined  and  the  average  value  of 
amplitude  and  phase  are  computed  for  the  included  points.  The  first 
subaperture  is  then  assigned  the  three  constants  Fx,  <;>i,  and  Lx  of 
Eq.  (17),  The  scan  continues  across  the  positive  half  of  the  aperture 
until  all  points  are  included,  returns  to  the  aperture  center  and  similarly 
scans  the  negative  portion  of  the  reconstructed  aperture.  Thus  the  n 
values  of  Fn(y),  4>n(y)  Ln  are  determined.  Equation  (17)  is  then 
evaluated  as  the  summation  of  N  integrals  having  uniform  illuminations. 
This  result  is  written  as: 


(17)  E(4.)  =  y 
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(18) 


E{<j>)  = 


F  n  Ln 


J4>n 


/kLn  sintj)'\ 

— j 

kLn  sintp 


where  the  term  <j)n  contains  an  additional  term  which  accounts  for  the  n-th 
subaperture  being  displaced  from  the  v'oordinate  axis.  Eulers  equation 
is  used  to  evaluate  Eq.  (18)  on  the  computer.  The  range  on  <{)  which  is 
calculated  is  pre-- as  signed  and  depends  upon  the  desired  end  result,  i.e., 
pattern  distortion  or  boresight  error.  The  method  of  scanning  the  aperture 
from  the  center  out  to  each  end  is  used  to  preseive  the  symmetry  of  the 
system. 


D.  Boresight  Error 

The  boresight  error  of  an  antenna- radome  system  can  be  defined 
as  the  difference  between  the  actual  target  direction  and  the  antenna 
pointing  direction.  In  a  well  designed  system  this  difference  is  a  few 
tenths  of  a  degree  and  is  due  primarily  to  phase  and  amplitude  distortions 
of  the  antenna  pattern  caused  by  the  radome.  The  boresight  error  is- 
evaluated  in  this  analysis  from  phase  monopulse  patterns  which  are 
generated  by  making  one-half  of  the  source  aperture  opposite  in  sign 
from  the  other  half.  This  pattern  is  characterized  by  a  deep  null  on  the 
beam  axis.  The  object  being  tracked  or  guided  by  the  particular  radar 
system  has  the  characteristic  direction  of  the  null  which  is  refered  to 
as  the  boresight  direction.  The  shift  in  the  location  of  this  null  due  to 
the  addition  of  a  radome  to  an  antenna  system  is  the  radome  boresight 
error.  If  the  antenna  is  scanning  in  a  partic-ular  direction  and  the  bore¬ 
sight  error  is  in  the  same  direction  it  is  defined  to  he  a  positive  error. 

In  calculating  the  boresight  error  several  considerations  simplify 
the  task.  The  null-shift  is  generally  a  fraction  of  a  degree,  thus  making 
the  calculation  of  only  a  small  portion  of  the  pattern  necessary.  Also, 
the  pattern  over  a  small  interval  enclosing  the  null  is  monotonically 
increasing  on  both  sides  of  the  null  and  approximately  symmetrical.  The 
null  location  is  determined  by  computing  one  pattern  point  on  each  side 
of  the  null  so  as  to  enclose  the  null  in  a  bracket.  By  use  of  the  symmetry 
and  monotone  properties  of  the  pattern  the  relative  values  of  the  two 
points  calculated  indicate  which  point  is  closest  to  the  null.  Fr^m  this 
information  a  third  point  is  calculated  which  halves  the  size  of  the  bracket 
containing  the  null.  Examination  of  the  field  magnitudes  at  each  end  of 
the  new  bracket  now  predicts  the  calculation  of  a  fourth  point  which  again 
halves  the  bracket  containing  the  null.  This  process  can  be  continued 
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indefinitely  to  obtain  the  null  location  to  any  desired  accuracy.  Starting 
with  a  two  degree  interval  the  null  location  will  be  known  to  within  1/2^ 
degrees  for  n  such  calculations.  In  this  analysis  an  n  of  11  is  used  which 
gives  an  accuracy  of  0.0005  degrees  in  the  null  location.  Figure  4 
illustrates  a  typical  calculation. 


NO  RADOME 
BORESIGHT  AXIS 


Fig.  4.  Far- zone  field  points  calculated  to  determine  the  null 
location  for  a  monopulse  difference  pattern  of  an 
antenna- radome  system.  The  order  of  the  points  cal¬ 
culated  3  indicated  by  the  number. 
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III.  DISCUSSION  OF  COMPUTER  PROGRAM 


The  computer  program  for  the  discussed  calculations  is  composed 
of  a  main  deck  and  several  subroutines  ?s  illustrated  in  Fig.  5.  The 
programs  are  written  in  Fortran  IV  language  suitable  for  calculation  on 
the  two  Ohio  State  University  computers,  the  IBM  7094  and  the  IBM 
System  360/75.  A  brief  description  of  the  function  of  each  routine  follows. 


MAIN  PROGRAM  SUBROUTINES 


F'g.  5.  Organization  of  computer  program. 
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1.  Main  Program:  Read  in  all  pertinent  data;  call  necessary  sub¬ 
routines;  calculate  antenna  patterns  with  and  without  radome; 
determines  relative  power  transmission  and  radome  bore- 
sight  error  as  a,  function  of  antenna  look  angle. 

2  .  Geometry:  determines  the  set  a,  b,  c,  d,  e,  f  of  geometry  coefficients 
for  each  of  the  n  radome  geometry  sections. 

3.  Ray  Trace:  determines  the  boundaries  between  the  n  radome  geo¬ 

metry  sections;  determines  the  n  X  m  matrix  of  incidence 
angles  corresponding  to  the  n  geometry  sections  and  the 
m  rays . 

4.  Multilayer  Transmission:  determines  the  transmission  coefficient 

and  insertion  phase  delay  for  the  n.  X  m  matrix  of  incidence 
angles . 

5.  Aperture  Taper:  determines  the  amplitude  and  phase  associated 

with  each  ray.  Also  calculates  any  aperture  blocking  or 
metal  nosecap  approximations. 

6.  Graph:  calculates  and  plots  the  normalized  far- zone  power 

pattern  in  dBs  with  and  without  radome. 

7.  Sidelobe  Level;  calculates  the  level  of  the  maximum  sidelobe  as 

a  percent  of  main  beam  intensity  and  as  dBs  down  from  main 
beam  intensity.  Also  calculates  the  location  of  the  first  side¬ 
lobe  for  the  no-radome  case  and  the  power  level  at  this 
location  with  the  radome  installed. 

8.  Half- Power  Beamwidth:  calculated  the  half-power  becimwidth  of 

the  antenna- radome  system  with  and  without  radome. 

Switching  from  the  main  program  to  the  desired  subroutines  is 
accomplished  by  means  of  two  input  cards  ncimed  "title"  and  "source" 
which  contain  key  words  describing  the  type  of  calculation  desired.  For 
example,  if  "no"  occurs  in  source  (3),  indicating  that  the  no-radome  case 
is  to  be  calculated,  the  multilayer  transmission  subroutine  -s  not  called. 
Comment  cards  have  been  placed  at  the  beginning  of  the  main  program 
which  explain  all  of  the  options  used.  Also,  most  program  variables 
are  explained  in  this  extensive  set  of  comment  statements. 

Figure  6  is  a  functional  flow  diagram  of  the  calculation.  The  signifi¬ 
cant  definitions  of  terms  used  in  this  diagram  are  listed  in  Table  I 
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Fig.  6.  Functional  flow  diagram  of  computer  calculation. 
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TAIfLE  I 


MC 

NOS 

N 

SHAPE 
(XOO,  YOO) 
ROO 
PHLIN 
(XLIN,  YLIN) 
AA,  ...FF 
LI. 

POLIZ 

A 

NUM 

DCE 

D 

TD 

FREQ 

PHD 

TRD 

SPAN 

PL 

RPL 

TAPER 

LT 

ANEA 

MI 

NRE 

NR 

AF 

CF 

TL,  CL,  SL 
CL2,SL2 
S2L 
MM 

XGB,  YGB 

YO 

LG 

TI 

NKRRI 

NKJRM 

AAT 


=  Number  of  cases  to  be  run. 

=  Number  of  geometry  sections 
=  Number  of  layers  in  each  section 
=  Geornetiical  shape  of  each  section 
=  Coordinates  of  center  of  an  ogive  section 
=  Radius  of  an  ogive  section 
=  Included  half-angle  of  a  cone  section 
=  Any  point  on  a  conical  section 
-  Geometry  constants  of  Eq.  (1) 

=  Total  number  of  lock  angles  used 
=  Polarization 
=  Length  of  source  aperture 

=  Number  of  points  calculated  in  partial  pattern 
=  Relative  dielectric  constant  of  a  layer 
=  Thickness  of  a  layer  in  inches 
=  Loss  tangent  of  a  layer 
=  Frequency  in  gigahertz 

=  Phase  allowance  uscd  in  numerical  integration 
=  Transmission  coefficient  allowance 
=  Angular  range  of  pattern  calculation 
=  J.ook  angle  in  degrees 
=  Look  angle  i  i  radians 
=  Definition  of  aperture  taper  used 
=  Present  value  of  LL 
=  Number  of  equal  length  subapertures 
=  ANEA 

=  Number  of  equal  spaced  rays 
=  NRE-1 

=  Fractional  length  of  a  subaperture  of  source 
=  Phase-center  correction  for  a  subaperture 
=  Tangent,  cosine,  sine  of  look  angle 
=  Cosine,  sine  of  twice  look  angle 
=  Sine- squared  of  Icok  angle 

=  Number  of  subaperture  mmediately  below  Y-axis 
(half- aperture  subdivision  point) 

=  Coordinates  of  geometry  bounds  between  geometry 
sections 

=  Y-coordiiiate  of  center  of  a  subaperture 
-•  Number  for  a  specific  geometry  section 
=  Angle  of  incidence 
=  Ray  trace  subroutine 
=  Multilayer  transmission  subroutine 
=  Aperture  amplitude  taper 
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TABLE  I  (Cent.) 


PHST 

= 

Aoerture  phase  taper 

FIPD 

= 

Insertion  phase  delay 

TC 

= 

Transmission  coefficient 

TR 

= 

Transmission  factor  for  a  subaperture 

PH 

= 

Phase  factor  for  a  subaperture 

RPWR 

= 

Relative  power  normalized  to  no  radorrie  case 

XMAX 

= 

Pattern  maximum 

BMAX 

= 

Angle  at  which  XMAX  occurs 

BNULL 

-2 

Angle  at  which  pattern  null  occurs 

BSEM 

= 

Boresight  error  in  milliradians 

SLL 

= 

Sidelobe  level  in  percent 

SLDB 

= 

Sidelobe  level  in  dB 

HPBW 

= 

Half-power  beamwidth  in  degrees 

To  indicate  the  execution  time  for  various  program  calculations 
the  following  table  is  presented.  500  rays  are  used  with  a  4  section 
radome  in  all  cases.  IBM  System  360/75. 


TABLE  II 


Calculation 


Execution  Time  for  500  Rays 


Ray  Trace 
Multilayer 
Taper 

Average  Aperture  Distribution 
Null 

100  Point  Partial  Pattern 
Plot  100  Pt  Pattern 
Combined  Sidelobe  and 

Half- Power  Beamwidth 


0.350  seconds 
0.6ii5  '  = 

0.083  ■' 

0.025  " 

0.050 

2.400  " 

0.250  ” 

0.017  ” 


It  is  seen  above  that  calculations  which  constitute  one  look  angle 
can  often  be  executed  in  less  than  one  minute. 
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IV. 


ANALYSIS  AND  DESIGN  OF 
ANTENNA- RADOME  SYSTEMS 


In  this  section  several  problems  in  antenna-iadome  system  design 
will  be  investigated  in  order  to  demonstrate  the  use  of  the  method  as  well 
as  to  point  out  its  applicability  to  a  wide  range  of  problems .  It  should  be 
emphasized  at  this  time  that  all  calculations  are  based  on  the  two-dimen¬ 
sional  model  of  the  antenna- radome  systorr:  and  that  the  accuracy  of  the 
calculations  is  unknown.  Verification  ot  results  is  possible  either  by 
comparison  with  measurements  or  by  comparison  with  results  obtained 
using  a  more  rigorous  theory.  A.s  was  stated  earlier  rigorous  theories 
presently  available  are  not  easily  applied  if  they  can  he  applied  at  all. 
Therefore,  whenever  possible,  results  will  be  compared  to  measured 
data. 


Two  specific  modern  radomes  configurations  will  be  used  in  most 
of  the  calculations  to  follow.  The  first  radome  is  characterized  as  a 
half-wave-wall  design  having  an  ogival  body  with  a  hemispheric  nosecap. 
The  aft  portion  of  the  radome  is  conically  faired  to  the  associated 
miissile  body.  Construction  is  entirelv  of  pyroceram  (cy  =  5.5),  The 
radome  wall  thickness  is  approximately  one-half  wavelength. 

The  fineness  ratio,  which  is  defined  as  the  ratio  of  the  axial  length  to 
the  base  diameter  of  a  radome,  is  2.0.  The  second  radome  is  derived 
from  the  first  by  removing  the  hemispheric  nosecap  and  extending  the 
ogive  body  to  form  a  closed  radome.  All  parameters  remain  the  same 
with  the  exception  of  the  Fineness  Ratio  which  becomes  2.25.  The  choice 
of  these  two  shapes  will  allow  an  evaluation  of  the  effects  of  blunting 
the  nose  of  a  ladcme,  which  is  frequently  necessary  because  of  aero¬ 
dynamic  heating  at  the  radome  tip.  Some  other  radom.e  configurations 
are  analyzed  which  will  be  specifically  described  v/hen  considered. 

Some  special  design  situations  require  modification  of  the  basic  method; 
these  will  be  pointed  out  when  necessary. 


A.  Convergence  Of  The  Ray- 
Optics  Solution 

As  w’as  pointed  out  in  Section  II- C  a  ray  tracing  calculation  has 
the  property  of  converging  to  a  fixed  answer  as  the  num.ber  of  rays  used 
is  increased.  This  section  presents  calculations  on  two  radome  geo¬ 
metries  to  illustrate  this  convergence  and  to  examine  the  number  of  rays 
required  to  obtain  the  convergent  solution.  Figure  7  shows  the  calcu¬ 
lated  boresight  error  versus  the  number  of  equally  spaced  rays  used  in 
the  calculation  for  the  pyroceram  radome  having  an  ogival  body  and  a 
hemispheric  nosecap.  Two  representative  look  angles  are  used  to 
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SOREStGHT  ERROR  (  MILLIRACIANS  ) 


I 


Fig*  ?-  Calculated  boreuight  error  vs.  number  of  equally  spaced 
rays  at  tVv’O  representative  look  angle.  Center  design 
frequency  aiid  perpendicular  polarization  for  the  pyroceram 
ogive  radojne  with  hemispheric  nosecap. 

illustrate  the  convergence;  complete  tabulatC'd  data  for  this  calculation 
at  ten  look  angles  are  included  as  Appendix  A..  Figure  8  shows  the  same 
calculation  for  the  radome  with  the  nosecap  removed  and  the  ogive 
extended  to  complete  the  radome.  Figure  9  shows  the  percent  difference 
from  the  final  answer  for  the  7.5*  look  angle  case.  It  is  seen  that  to 
obtain  the  convergent  solution  {Ofc  error)  a  large  numbex'  ot  equally- 
spaced  rays  is  required.  The  presence  of  the  nosecap  is  seen  to  have 
little  effect  on  the  convergence  of  the  solution  if  more  than  10  rays  are. 
used. 
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Fig.  8.  Calculated  boresight  error  vs.  number  of  equally  spaced 
rays  at  two  representative  look  angles.  Center  design 
frequency  and  perpendicular  polarization  for  the  pyro- 
ceram  ogive  radome  without  hemispheric  nosecap. 


B .  Numerical  Integration  Of  The 
Reconstructed  Aperture 

This  section  demonstrates  the  convergence  obtained  using  the 
numerical  integration  technique  of  weighted  subaperbares  explained  in 
Section  II-C.  501  equally- spaced  rays  are  used  to  represent  the  ten 
wavelength  aperture  used  with  the  radomes  of  Figs.  7  through  9  in  all 
of  the  following  calculations.  Boresight  error  (BSE),  relative  on-axis 
power  (RPWR),  and  the  number  of  weighted  subapertures  obtained  (N) 
are  calculated  for  various  rombinations  of  phase  allowance  (PHD)  and 
amplitude  allowance  (TRD)  in  approximating  the  field  over  each  sub¬ 
aperture  by  uniform  amplitude  and  phase.  Table  III  shows  some  lepre- 
sentative  calculated  results  with  BSE  and  PHD  in  degrees.  Appendix  B 
gives  the  complete  data  for  this  calculation  at  ten  look  angles.  The 
values  obtained  for  0  amplitude  and  0®  phase  allowances  are  the  same 
results  obtained  in  Section  A  of  this  Chapter,  i.e.,  the  convergent  answer 
from  a  large  number  of  equally  spaced  rays.  The  nuinber  of  aperture 
points  indicated  in  the  table  are  the  number  of  subsections  of  the  aperture 
which  result  for  a  given  phase  and  amplitude  allowance  combination  and 
indicate  the  relative  time  consumption  for  a  computer  pattern  calculation. 
The  table  shows  that  the  convergent  answer  is  obtained  using  almost  any 
of  the  given  allowances  -  for  the  complete  ogivo  radome  the  answer  is 
obtained  using  as  few  as  7  subsections  of  the  original  500  point  aperture 
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E 


Fig.  9.  Percent  difference  from  convergent  ray-tracing  solution 
vs.  nximber  of  rays.  Perpendicular  polarization  at 
center  design  frequency.  Look  angle 


approximation.  With  the  large  phase  and  amplitude  allowances,  30“  and 
0.2  respectively,  the  error  is  still  less  than  2%.  In  calculations  to  follow 
an  amplitude  allowance  of  0.1  and  a  phase  allowance  of  3 . 0  degrees  are 
generally  used.  501  e4usr>  spaced  rays  are  used  throughout  since  this 
number  assures  that  the  vonvergent  solution  can  be  obtained.  The  con>- 
bined  time  consumption  for  tracing  501  rays  and  computing  the  associated 
values  of  jx!*  and  IPD  is  only  on  the  order  of  1  second.  As  was  pointed 
out  earlier  the  significant  computer  time  usage  occurs  in  calculating 
pattern  points.  This  is  because  each  pattern  point  requires  summing 
contributions  from  each  aperture  point  used.  For  example  using  500 
equally  spaced  rays,  10  look  a.ngles,  and  calculating  100  points  in  the 
far-field  pattern  the  total  number  of  calculations  is  500,  000.  If  only  10 
aperture  points  are  used  the  munber  of  calculations  is  reduced  to  10,  000 
which  is  quite  significant. 
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C.  Electrical  Design  Of  A  Radome  Wall 


Tlie  high  speed  attained  using  the  two-dimenoicnal  analysis  results 
in  relatively  low  cost  calculations.  This  allows  the  method  to  be  used 
to  advantage  as  a  design  tool.  The  approach  is  to  select  cin  approximate 
design  in  terms  of  the  complex  dielectric  constants,  wall  thicknesses, 
number  of  layers,  and  geometrical  shape.  A  specific  parameter,  for 
example  wall  thickness,  is  varied  in  small  steps  above  and  below  the 
design  specification.  Calculations  of  desired  electrical  parameters,  such 
as  boresight  error,  transmission,  sidelobe  level,  etc.  are  made  at  each 
incremental  variatic  t.  Data  are  then  compared  to  determine  an  optimum 
design.  One  such  excimple  follows. 

Figure  10  shows  the  calculated  boresight  error  for  the  pyroceram 
radome  as  a  function  of  its  wall  thickness.  The  current  design  thick¬ 
ness  is  specified  as  0  percent.  From  the  curves  it  can  be  seen  that 
for  any  look  angle  the  boresight  error  approaches  a  low  value  in  the 
range  of  -3  to  -5  percent.  Further,  in  this  range  the  actual  value  of 
error  for  a  given  look  cuigle  remains  r  atively  constant.  This  indicates 
that  the  radome  would  operate  well  in  uiis  region  and  show  practically  no 
change  in  error  due  to  small  frequency  drifts,  dimensional  tolerances,  or 
thermal  gradients . 

If  we  examine  the  curves  near  0  percent  or  higher  v/e  find  that  the 
radome  will  be  sensitive  to  the  above  three  mentioned  considerations 
and  operate  with  significantly  higher  boresight  error  as  well.  Thus  it 
appears  that  a  4  percent  reduction  in  wall  thickness  would  reduce  the 
maximum  boresight  error  by  50  percent.  The  on- axis  transmission 
efficiency  of  the  radome  is  calculated  simultaneously  with  the  boresight 
error  in  order  that  the  effects  of  a  design  change  on  transmission  can  be 
observed.  Figure  11  show's  that  decreasing  the  wall  thickness  by  4  per¬ 
cent  causes  a  14  percent  net  loss  in  transmitter  power.  This  is  probably 
not  excessive  in  view  of  the  improvement  in  boresight  performance. 


D.  Radome  Boresight  Error  Versus 
System  Bandwidth 


If  the  calculated  curves  of  Fig.  10  are  correct,  precise  agreement 
between  calculated  and  measured  da.ta  in  the  region  of  design  thickness 
(0%)  ife  unexpected.  A  small  dimensional  error  could  easily  cause  a  15- 
20  percent  change  in  boresight  error.  Figure  12  shows  a  comparison 
between  calculated  and  measured  data  for  the  pyroceram  radome.  Agree¬ 
ment  is  only  fair  in  this  case.  The  frequencies  in  Fig.  12  correspond  to 
the  upper  and  lower  frequencies  of  a  1.5  percent  bandwidth  design. 
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Fig.  10.  Calculated  boresight  error  for  pyrocerarn  ogive 
radome  with  hemispheric  nosecap  as  a  function 
of  its  wall  thickness.  Perpendicular  polariza¬ 
tion  at  design  frequency. 


E .  Source  Taper  Effects 


Many  radar  designs  utilize  carefully  controlled  amplitude  tapers 
in  order  to  achieve  an  antenna  pattern  having  very  low  sidclobes.  Phased 
array  techniques  available  today  emphasize  this  method.  A  studv  to 
determine  the  pattern  distortion  in  terms  of  change  in  sidelobc  level  and 
half-power  beamwidth  due  to  the  addition  of  a  radome  to  such  an  antenna 
system  was  carried  out.  In  addition,  the  effects  of  the  use  of  an  ampli¬ 
tude  taper  on  the  system  boresight  error  characteristics  were  calculated. 
Two  antenna- radome  systems  were  analyzed,  the  pyrocerarn  ogive  radome 
with  the  hemispheric  nosecap  {blunted  nose  case)  and  the  complete  ogive 
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Fig.  1)  .  On-axis  transmission  efficiency  for  pyroceram  ogive 
radome  with  hemispheric  nosecap  as  a  function  of  its 
wall  thickness.  Perpendicular  polariz.ation  at  design 
frequency. 


radome  (pointed  nose  case).  An  identical  antenna  having  a  variable 
amplitude  taper  was  analyzed  for  the  two  radomes.  Particular  emphasis 
was  placed  on  the  "cosine-squared  on  a  pedestal"  distribution  since  it 
provides  a  convenient  method  for  varying  the  antenna  pattern  over  a 
broad  range  of  sidelobe  levels.  Also,  this  distribution  is  commonly 
used  to  achieve  low- sidelobe  pencil  beam  antennas. 

Figure  13  shows  the  sidelobe  level  obtained  using  the  two  antenna- 
radome  systems.  The  blunted  nose  radome  is  seen  to  degrade  the  antenna 
performance  severely  while  the  pointed  nose  radome  produces  inconse¬ 
quential  pattern  distortion.  All  calculations  were  made  at  look  angle  0“ 
in  order  to  emphasize  the  difference  between  the  two  systems.  These 
calculations  are  corroborated  to  som.e  extent  by  measurements  performed 
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BORESIGHT  ERROR  (MILLIRADIANS) 


Fig.  12.  Calculated  and  measured  boresight  error  for  pyroceram 
ogive  radome  with  hemispheric  nosecap.  F’jigh 
Flov/  denote  the  two  extremes  of  a  1.5%  bandwidth 
design.  Perpendicular  polarization. 


on  a  similar  antenna -radome  system  by  Styron  and  Hoots®  of  the 
Brunswick  Corporation.  They  measured  pattern  distortion  due  to  a 
blunted  nose  conical  radome  in  terms  of  sidelobe  degradation  for  three 
aperture  tapers.  They  found  that  a  basic  30  dB  sidelobe  antenna  was 
reduced  to  an  approximately  21  dB  system  and  that  the  radome  controlled 
the  sidelobe  level  rather  than  the  aperture  taper.  No  similar  set  of 
measurements  are  available  for  the  pointed  nose  radome,  however, 
Styron  a.id  Hoots  stated  that  the  pattern  degradation  was  most  severe 
at  offsets  where  the  center  antenna  ray  impinged  near  the  radome  nose 
and  that  for  offsets  further  from  the  nose  the  degradation  was  minimal. 
This  tends  to  verify  the  calculations  for  the  pointed  nose  case. 
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SIDELOBE  LEVEL  (dB) 


PEDESTAL  HEIGHT 


Fig.  13.  Calculated  sidelobe  level  for  the  pyroceram  ogive 
radon^e  with  and  without  hemispheric  nosecap. 
Look  angle  0®,  perpendicular  polarization,  design 
frequency.  Aperture  amplitude  taper  is  cosine - 
squared  on  a  variable  pedestal. 
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Figure  14  shows  that  the  half-power  beamwidth  is  relatively  un¬ 
affected  by  the  presence  of  the  blunt  nose  radome.  Similar  results 
were  obtained  for  the  pointed  nose  case. 

Figures  15  and  16  shov/  the  effects  of  several  amplitude  tapers  on 
the  boresight  error  performance  of  the  two  antenna- radome  systems^ 
The  blunted-nose  radome  is  seen  to  have  considerably  poorer  boresight 
performance  when  an  amplitude  taper  is  used.  This  may  be  attributed 
to  the  much  smaller  radius  of  curvature  in  the  vertex  region  which 
causes  considerable  phase  distortion  in  the  aperture  distribution.  Also, 
since  this  ray-tracing  analysis  uses  a  collimated  beam  projecting  from 
the  source  through  the  radome,  it  is  likely  that  the  resulting  higher  con¬ 
centrations  of  enerey  near  the  vertex  region  tend  to  over-emphasice  the 
effects  of  the  ver:ex  region.  Thus  the  effects  of  the  blunted  nose  on  the 


Fig.  14.  Effect  of  an  aperture  amplitude  taper  on  antenna- 
radome  system  bandwidth.  Blunt  nose  pyroceram 
ogive  radome,  design  frequency,  perpendicular 
polarization. 
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Fig.  15.  Effects  cf  several  aperture  amplitude  tapers  on 
system  boresight  error.  Pointed  nose  pyro- 
ceram  ogive  radOiiie,  perpendicular 
polarization  at  design  frequency. 


radome  are  somewhat  exaggerated.  From  the  error  curves  it  is  seen 
that  in  both  cases  the  more  nearly  uniform  aperture  distributions  pro¬ 
duce  the  lowest  boresight  error.  The  best  boiesight  performance  ’s 
obtained  with  the  uniform  distribution,  however,  the  cosine- squared  on 
a  pedestal,  or  something  similar,  can  be  used  for  sidelobe  control  with- 
o\it  seriously  affecting  the  boresight  performance.  Measurements  by 
Styron  and  Hoots*  support  the  above  calculations  for  the  blunted  nose 
radome  case. 


F .  Aperture  Blocking 

In  the  two  dimensional  ray-tracing  approximation  a  metallic  portion 
of  a  radome,  such  as  a  protective  rain-erosion  cap,  is  treated  as  an 
aperture  block.  This  requires  specific  changes  in  the  computer  program 
for  two  reasons.  First,  the  perfectly  collimated  beam  assumed  in  the  ray 
tracing  approach  predicts  that  the  effects  of  ar.  obstacle  in  front  of  an 
antenna  are  independent  of  the  distance  between  the  obstacle  and  the 
antenna.  The  second  problem  is  that  the  portion  of  the  source  aperture 
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BORESIGHT  ERROR  (  M  ILLIRADI ANS) 


LOOK  AMGLE  (  DEGREES  ) 


Fig.  16.  Elfect-D  of  several  aperture  amplitude  tapers  on 

system  boresight  error.  Pointed  nose  pyroceram 
ogive  radonie,  per pendicular  polar- nation  at 
design  frequency. 
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blocked  in  the  two-dimensional  model  is  m  ich  larger  than  the  actual 
area  blockage  in  ‘■he  thrv-e-dimensional  problem.  A  study  was  made  to 
determine  a  suitable  two-dimensional  representation  of  the  three- 
dimensional  block.  Details  of  the  specific  treatment  for  aperture 
blocks  are  given  below. 

Figure  17  shows  an  aperture  block  of  radius  h  located  at  a  distance 
/  from  an  antenna  aperture  of  radius  R.  The  source  aperture  is  pro¬ 
jected  to  the  plane  of  the  block  using  a  divergence  angle  Gj  equal  tc  the 
half-power  bcamwidth  to  determine  a  projected  aperture  radius  R': 


(19)  R'  =  R  +  £  sin  Qi 


The  ratio  of  the  block  area  to  the  projected  source  area  is  calculated 
as; 


Fig.  17.  Geometry  ■*  aperture  blocking  calculation. 
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(20)  Ratio  =•  - - - 

(R  +  f  sinOi)^ 

Equation  (20)  gives  that  fraction  of  the  source  area  to  be  blocked  out 
for  any  value  of  £  .  In  the  one-dimensional  aperture  approximation  for 
a  planar  source  the  block  is  inserted  at  £  =  0  even  though  the  block  is 
located  at  t ;  hence  h  must  be  reduced  to  account  for  the  distance  £  . 
The  effective  blocked  area  at  the  source  is: 


(21)  -Ab LOCKED  =  -AsOURCE  * 


which  gives: 


(22) 


h'  = 


BLOCKED 

ir 


) 


2 


as  the  reduced  length  of  the  block.  This  is  the  approximate  method  used 
to  account  for  the  antenna- obstacle  separation. 

A  second  approximation  is  required  because  the  one-dimensional 
block  does  not  represent  the  two-dimensional  block  in  the  other  uimension. 
Figure  18  shows  that  the  blockage  in  the  two-dimensional  case  represented 
by  the  one-dimensional  block  is  a  strip  across  the  entire  aperture.  The 
approximation  used  here  is  to  reduce  the  block  length  such  that  the  re¬ 
sulting  strip  area  is  equal  to  the  actual  area  of  the  block.  In  this  way# 
even  though  the  shapes  of  the  aperture  blocks  differ,  the  source  area 
blocked  out  is  the  same.  With  reference  to  Fig.  19,  the  str’p  area  is: 


(23)  AgTRIP  =  R^(tt  -  6  +  sin0) 


where  0  in  radians  is  given  by: 


(24)  0  =  2  cos”'  (y/R) 


Using  Eq.  (21)  we  set 
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Fig.  18.  Effective  aperture  blockage  in  two  dimensions  by 
a  one-dimensional  aperture  block. 


(25)  ttR*  X  RATIO  =  R*(tt  -  0  -  sin 9) 

Removing  the  R^  terms  and  rew'riting  Eq.  (25)  in  homogeneous  form; 

(26)  9-sinG  ir  ( 1  -  RATIO)  =  0 


This  equation  can  be  solved  to  any  degree  of  accuracy  using  Newton's 
method 


(27) 


f(a,) 

a?  =3' - : — i — 

‘  f’(a,) 
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Fig.  19 «  Geometry  used  in  calculating  the  area  of 
a  circular  strip. 


where  ai  is  an  approximate  solution  of  Eq.  (26)  which  is  f  in  Eq.  (27), 
The  value  3.z  is  a  better  approximate  solution  than  aj.  By  iterating 
Eq.  (27)  we  can  obtain  any  desiied  accuracy  in  the  approximation. 
Since  most  blocks  are  small  a  value  of  170®  is  used  for  aj  in  the  pro- 
“gram.  In  the  problem  being  considered  Eq.  (27)  takes  the  form 


e,  =  8, 


Gj  -  sinOi-  Tr(l  -  RATIO) 
1  -  cos  Oj 


The  resulting  block  width  in  the  one-dimensional  aperture  approximation 
using  this  approximation  is:' 


2  R  cos(0/ 2) 


where  0  is  the  angle  associated  with  the  strip  as  shown  in  Fig.  19. 

In  order  to  determine  the  accuracy  obtained  using  the  above  two- 
dimensional  aperture  blocking  approximations  some  Scimple  calcula¬ 
tions  were  made  and  compared  to  the  calculated  and  measured  aperture 
blocking  results  of  Collier.*  In  his  report  Collier  used  a  16.4X.  para¬ 
bolic  dish  having  a  2.2X.  diameter  feed  located  2.2X  in  front  of  the 
aperture  plane.  He  considered  obstacles  ranging  in  size  from  3,4\ 
to  9.7X,  ^vhich  could  be  positioned  from  20\  to  lOOX.  from  the  aperture 
plcure.  The  frequency  v/as  32.7  GHz.  The  medium  sized  obstacle  of 
4.7X.  was  chosen  for  comparison  here.  Collier  used  a  severe  radial 
>taper»  as  shown  in  Fig.  20,  to  represent  the  antenna  aperture  distri¬ 
bution.  The  central  amplitude  of  zero  was  used  to  account  for  the 
aperture  blocking  due  to  the  feed.  The  measured  pattern  reported 
showed  an  approximately  5*  beamwidth  and  16., 5  dB  sidelobes.  This 
taper  was  represented  here  by  a  piece'' ise  linear  approximation  with 
the  exception  that  the  curve  was  extended  to  1  0  at  the  origin  and  the 
feed  treated  as  a  separate  aperture  block.  The  calculated  pattern  using 


Fig.  20.  Aperture  amplitude  taper  used  in  Collier  aperture 
blocking  calculation. 
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the  two-dimensional  approximations  showed  a  beamwidth  of  5.29*  and  a 
sidelobe  level  of  16.32  dB  which  is  in  excellent  agreement  with  Collier's 
measurements.  Figure  21  shows  the  calculated  sidelobe  level  for  the 
antenna  in  the  presence  of  the  obstacle  as  a  function  of  aperture-obstacle 
separation.  The  modified  two-dimensional  model  generally  shows  very 
good  agreement  with  measurements  and  with  the  three-dimensional 
calculations . 


C- .  Electrical  Performance  Of  A  Radome 
In  A  Hyper -Environment 

Due  to  the  high  speed  of  modern  aircraft  and  missiles,  radomes 
are  often  subjected  to  severe  environments.  Nonuniform  temperature 
distributions  exist  about  the  radome  wall  which  result  in  variations  in 
the  temperature  dependent  quantities  of  dielectric  constant,  loss  tan¬ 
gent,  and  wall  thickness.  The  variations  in  these  quantities  alter  the 


bovesight  error  performance  of  the  antenna- radome  system.  Figure  22 
shows  a  representative  temperature  profile  which  a  radome  may  en¬ 
counter.  The  boresight  error  characteristics  of  the  pyrccsram  radome 
having  a  hemispheric  nosecap  were  calculated  for  this  temperature 
profile  as  an  example.  To  approximate  the  effects  of  the  temperature 
profile,  the  radome  is  subdivided  into  several  sections,  each  of  which 
has  a  fixed  set  of  dielectric  constant,  loss  tangent,  and  wall  thickness 
parameters.  The  subdivision  is  determined  by  observing  tha  rate  of 
change  of  these  parameters  with  temperature  and  the  rate  of  change  of 
temperature  along  the  radome  wall.  References  7  and  8  were  used 
for  this  purpose.  Figure  23  shows  the  calculated  boresight  error  in 
the  presence  of  the  temperature  profile  of  Fig.  22.  Figure  24  shows 
a  similar  set  of  curves  with  the  original  wall  thickness  reduced  by 
1.5%.  These  results  indicate  that  the  boresight  error  performance  of 
the  radome  may  actually  improve  in  a  severe  thermal  environment. 

The  effect  of  the  temperature  profile  in  this  case  is  seen  to  be  similar 
to  the  design  technique  of  constructing  a  tapered  radon,?,  wall  to  improve 
radome  performance. 

As  shown  in  Fig.  22  there  is  a  temperature  gradient  from  the 
outside  to  the  inside  of  the  radome  wall.  In  the  above  example  the  tem¬ 
perature  was  assumed  to  vary  linearly  with  distance  through  the  wall. 

In  case  there  is  a  nonlinear  variation  in  temperature  or  parameter 
constants  as  a  function  of  temperature  through  the  wall,  a  fuither  approxi¬ 
mation  consisting  of  subdividing  each  section  into  several  layers  having 
variable  parameters  can  be  used.  Thus  the  final  subdivision  of  the  radome 
in  this  case  would  be  one  of  several  geometry  sections  having  differing 
numbers  of  layers . 


H.  Comparison  Of  Boresight  Measurements 
And  Calculations 


Boresight  calculations  using  the  ray  tracing  method  described  in 
this  report  are  compared  in  this  section  with  measured  data  supplied  by 
The  U.S.  Naval  Air  Development  Center  at  Johnsville,  Pa.  and  with 
calculations  and  measurenients  taken  from  the  literature.  Exact 
radome  geometry  was  not  always  known  in  the  following  cases,  hence 
some  comparisons  were  made  on  the  best  estimate  basis. 

Case  1 


The  previously  described  half-wave  wall  conically-faired  ogive 
radome  with  a  hemispheric  nosecap  is  examined  here.  A  constant  wall 
thickness  of  pyroceram  (e  =  5.5)  was  used  throughout.  Calculations  and 
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Fig.  23.  Calculated  boresight  error  for  bluni--nosed  pyroceram 
ogive  radome  in  Mie  presence  of  the  temperature  pro¬ 
file  of  Fig.  22.  Perpendicular  polarization  at  desigh 
frequency. 
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Fig.  24,  Calculated  boresight  error.  Same  case  as  Fig.  23 
with  radome  wall  thickness  reduced  by  1.5%. 
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measurements  at  the  high,  low,  and  center  frequencies  of  a  1.5%  bai  d- 
width  design  were  made  for  both  perpendicular  and  parallel  polari-  ^ 
zations.  Two  sets  of  measuremern.s  were  furnished  by  the  USNADC;  ’ 
botli  cf  which  are  included  in  the  comoarvsons  to  give  an  indication  o: 
experimental  deviations.  This  deviation  is  generally  a  result  of  a 
lack  of  symmetry  in  the  radome.  Figures  ?5  and  26  show  the  compari¬ 
sons  between  calculations  and  measurements  for  perpendicular  and 
parallel  polarization  at  the  low  Irequcncy  end  of  the  band.  Agreement 
between  calculation  and  measurement  s  reasonably  good  for  the  per¬ 
pendicular  polarization  case,  however,  the  agreement  is  poor  for 
parallel  polarization.  Figures  27  and  28  give  the  same  comparisons 
for  the  design  frequency.  Perpendicular  polarizaticn  shows  very  gooo 
agreement  in  this  case  while  agreement  remains  poor  for  parallel 
polarization.  Figures  29  and  30  show  the  comparison  for  the  high  end 
of  the  frequency  band.  Agreement  is  also  very  good  for  the  case  of 
perpendicular  polarization  and  poor  for  parallel  polarization. 


LOOK  ANGLE  (DEGREES) 


Fig.  25.  Calculated  and  messured  boresight  error  for  a 
blunted-nose  pyrocereim  ogive  radome  at  the 
lower  limit  of  a  1.5%  bardwidth  design.  Per¬ 
pendicular  polarization. 
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Fig.  26.  (Tai._ulated  and  measured  boresight  error  for  a 
blunted-nose  pyroceram  ogive  radome  at  the 
lower  limit  of  a  1.5%  bandwidth  design.  Parallel 
polarization. 

Case  2 

This  comparison  uses  the  same  radome  geometry  as  in  Care  1 
except  that  the  construction  is  of  polyimide  («  =4.2).  Measurements 
and  calculations  at  the  design  frequency  for  perpendicular  and  parallel 
polarizations  are  shown  in  Fig.  3l.  Agreement  is  good  for  both 
polarizations  in  this  example.  Measurements  were  furnished  by 
USNADC.'° 
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Fig.  27.  Calculated  and  measured  boresight  error  for  a 
blunted-nose  pyroceram  ogive  radome  at  the 
design  frequency.  Perpendicular  polarization. 


Case  3 


This  comparison  uses  a  radome  having  basically  the  same  con¬ 
struction  as  in  Case  1,  i.e.,  conical  fairing,  ogive  body,  and  hemi¬ 
spheric  rosecap.  The  main  body  of  the  radome  was  constructed  of 
polyimide  (e  =  4.2).  A  rain  erosion  cap  of  alumina  (c  =  8.9)  was  sprayed 
on  the  tip  end  of  the  radome  extending  back  six  inches.  Exact  dimensional 
data  was  unavailable  for  this  case,  therefore,  estimates  were  made  as  to 
the  probable  design.  Calculations  of  boresight  error  (BSE)  and  tran.s- 
mission  efficiency  (It|^)  were  m.ade  at  the  probable  design  thickness  and 
at  i  2.57o  and  +  5%  increments  about  this  design  thickness.  In  the  region 
of  the  alumina  nosecap  the  main  body  of  the  radome  was  assumed  to  be  of 
thinner  wall  to  adjust  for  the  increased  dielectric  constant  of  the  cap 
material.  The  source  of  measured  data  was  a  report  by  The  Brunswick 
Corporation. “  Reference  11  also  stated  that  the  radome  had  been  "cor¬ 
rected"  -  which  amounts  to  adding  patches  of  dielectric  material  to  the 
radome  inner  wall  to  locally  alter  the  i;hase  shift  value  through  the  wall. 
Thus  the  geometry  is  relatively  uncertain.  High  and  low  frequency 
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Fig.  28.  CalculPted  and  measured  boresight  error  for  a 
blunted-nose  pyroceram  ogive  radome  at  the 
design  frequency.  Parallel  polarization. 


, calculations  were  made  at  perpendicular  polarization  since  these  corre¬ 
sponded  to  the  reported  measurements.  Figures  32  and  33  show  the  high 
and  low  frequency  boresight  errer  calculations  for  the  estimated  geo¬ 
metrical  construction  compared  to  the  measurements  reported  in  Ref.  11. 
Measurements  were  made  at  several  roll  angles  resulting  in  the  spread 
of  value  shown  by  the  shaded  portions  in  Figs.  32  and  33.  Agreement 
was  not  expected  to  be  good  in  this  case,  however,  consideruig  the 
assumptions  required,  agreement  is  satisfactory.  The  general  shape 
of  the  curves,  i.e.,  the  initial  negative  error  at  small  look  angles  and 
the  shift  to  positive  error  at  larger  look  angles  is  predicted.  The  ampli¬ 
tude  of  the  negative  swing  is  considerably  larger  than  the  measurement, 
however,  the  positive  error  an-.pi-’tude  is  in  good  agreement  with  the 
measurements.  This  could  be  explained  by  the  ''correction"  performed 
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Fig.  29.  Calculated  and  measured  boresight  error  for  a 
blunted-noss  pyroceram  ogive  radome  at  the 
upper  limit  of  a  1.5%  bandwidth  design.  Per¬ 
pendicular  polarization. 


on  the  radome  as  mentioned  above.  Using  the  main  body  (aft  of  the 
noGecap)  wall  dimensions  corresponding  to  Figs.  32  and  33  as  a  base 
the  radome  error  was  then  calculated  at  +  2.5%  and  ±  5%  of  this  value 
with  the  nose  dimensions  held  constant.  In  all  cases  the  rain  erosion 
cap  was  taken  to  be  0.030  inches  thick.  Figures  34  und  35  shew  these 
calculations  compared  to  the  measurements.  Excellent  agreement 
between  calcvilatioris  and  measurements  is  obtained  fc'  the  +  2.5%  case 
incidating  that  this  wall  thickness  is  probably  closer  to  the  actual  value 
than  the  original  assumed  value.  Another  view  of  this  same  data  is 
shown  in  Fig.  37  which  indicates  that  +  2.5%  is  about  optimum  for  a  bore 
sight  error  design  point.  The  computer  program  simultaneously  calcu¬ 
lates  transmission  efficiency  with  boresight  error.  Figure  37  shows  the 
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Fig.  30.  Cr.lculaci’d  and  measured  boresight  error  for  a 
blunted-nose  pyroceram  radome  at  the  upper 
limit  of  a  1.5%  bandwidth  design.  Parallel 
polar-!  zaiion. 

It]*  curves  which  correspond  to  Figs.  33  through  36.  Since  radomes 
cf  this  type  are  frequently  designed  for  maximum  transmission,  the 
indication  is  that  the  design  wall  thickness  was  +  2.5%  or  more  above 
the  assumed  value. 


Case  4 

Several  radome  body  geometries  previously  analyzed  by  General 
li)ynaniics  are  examined.  Calculations  of  boresight  eri  r  are  rr.'ade 
and  compared  to  similar  calculations  by  G.P.  Tricoles  of  General 
Dynamics  Corporation,  San  Diego,  California.*^  Since  Tricoles  uses 
a  different  ray-tracing  calculation  these  ccmparisons  are  presented 
merely  to  show  similarities  and  disagreements  between  the  calculated 
results.  Table  IV  indicates  the  geometrical  differences  among  the 


Fig.  31.  Calculated  and  measured  boresight  error  for  a 
blunt-nose  polyimide  radome  at  the  design 
frequency. 


various  radome  configurations.  Table  V  follows  the  radorne  curves  to 
indicate  I'he  general  agreement  between  the  two  calculations.  The  com¬ 
parisons  are  given  for  the  several  configurations  at  the  h'gh  and  low 
ends  of  a  1.5%  bandwidth  for  perpendicular  and  parallel  polarizations 
in  Figs.  38  through  51  . 
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Fig.  34.  Comparison  batween  measurements  and  calculations 
at  i  2.5%  and  t:5%  of  the  probable  design  thickness 
of  the  main  body  .or  the  ra.dome  of  Fig.  32.  Per¬ 
pendicular  polarization,  high  frequency. 
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Fig.  35.  Comparison  between  measurements  and  calculations 
at  +2.  5%  and  ^5%  of  the  probable  design  thickness  of 
the  main  radome  body  for  the  radome  of  Fig.  32. 
Perpendicular  polarization  at  the  low  fiequency. 
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Fig.  37.  Calculated  transmission  efficiency  vs  wall 
thickness  vs  look  angle  for  the  radome  of 
Fig.  32,  Perpendicular  polarization  at  the 
high  frequency. 
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TABLE  IV 

Description  of  Seven  Example  Radomec,  All  Dimensions 
are  in  Inches.  Sketch  Pertains  to  Secant  Ogives  Only. 


No. 

Fine¬ 

ness 

Ratio 

Shape 

1 

1 

2.0 

Tangent  Ogive 

5.5 

2 

2.11 

Tangent  Opive 

5.5 

3 

2.0 

Tangent  Ogive 

6.4 

4 

2.0 

Tangent  Ogive 

9.7 

5 

2.0 

Secant  Ogive 

5.5 

6 

2.5 

Secant  Ogive 

5.5 

7 

2.5 

Secant  Ocive 

5.5 

Geometrical  Constants 


R 


-1 

-1 

-1 


30.0 

37,5 


27.5 

95.0 

95.0 


5.692 

9.35 

9.35 


37.3 


119.031 

186.529 


KILCOYNE  CALCULATION  FH 
KILCOYNE  CALCULATfON  FL 


Comparison  between  Tricoles  and  Fig.  39.  Comparison  between  Tricoles  and 

Kilcoyne  calcv\lat’ons .  Radome  1  Kilcoyne  calculations.  Radome  1 

in  Table  IV.  Parallel  polarization.  in  Table  IV.  Perpendicular  polari 

FXj  and  Fh  denote  the  low  ai.d  high  zation. 

freouencies  of  a  1.5%  bandwidth. 


Comparison  between  Tricoles  and  Fig.  41.  Comparison  between  Tricoles  and 

Kilcoyne  calculations.  Radome  2  Kilcoyne  calculations.  Radome  2 

in  Table  IV.  Paraillel  polarization.  in  Table  IV.  Perpendicular 

polarization. 


XKILCOVNE  CAt-CULAFION  FH 
aKILCOYNE  CALCULATtON  FL 


iKILCOYNE  CALCULATION  FH 
KILCOYNE  CALCULATION  FL 
TRICOLES  CALCULATION  FH 
TRICOLES  CALCULATION  FL 


Comparison  between  Tricoles  and  Fig«  45.  Comparison  between  Tricoles  and 

Kilcoyne  calculations.  Radome  4  Kilcoyne  calculations.  Radome  4 

in  Table  IV.  Parallel  polarization.  in  Table  IV.  Perpendicular 
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Comparos in  between  Triroles  and  Fig.  47.  Connparis on  between  Tricoles  a 

Kilcoy  "ve  calculations  .  Radomr  5  Kile oyne  calculations .  Radome 

in  Table  IV.  Parallel  polarization.  in  Table  IV.  Perpendicular 
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COYNE  CALCULATION 
COYNE  CALCULATION 
ICOLES  CALCULATIO^ 
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Fig.  50.  Comparison  between  Tricoles  and  ‘  Fig.  51.  Comparison  between  Tricoles  and 

Kilcovne  calculations.  Radome  7  Kilcoyne  calculations.  Radome  7 


TABLE  V 

Relative  Agreement  Bet'.veen  Tri-oles  and  Kilcoyne  Calculations 

FH  =  1.015  FL 


Radome 

Perpendicular 

Parallel 

FR 

Shape 

€ 

FL 

FH 

FL 

FH 

2.0 

Tan- Ogive 

5.5 

Poor 

Fair 

Poor 

Poor 

2.11 

Tan- Ogive 

5.5 

Good 

Excellent 

Excellent 

Excellent 

2.0 

Tan- Ogive 

6.4 

Poor 

Poor 

Fair 

Poor 

2.0 

Tan- Ogive 

9.7 

Poor 

Fair 

Good 

Poor 

2.0 

Sec  -Ogive 

5,5 

Poor 

Excellent 

Poor 

Poor 

2.5 

Sec  -  Ogive 

5.5 

Poor 

Fair 

Good 

Poor 

2.5 

Sec  -  Ogive 

5.5 

Poor 

Excellent 

Excellent 

Poor 

Weighted  Average 

Poor 

Good 

Good 

poor 

I.  Effects  IXie  to  the  Blunting  of  a 
Radome  Nose  Section 


The  high  speed  of  a  modern  aircraft  or  missile  frequently  results 
in  the  generation  of  temperatures  at  the  leading  edge  (tip)  that  are  above 
the  maximum  safe  operating  temperatures  of  even  the  best  ceramic 
radome  materials.  When  it  is  necessary  to  locate  a  radome  at  the  tip 
sectioni  special  design  precautions  are  required  which  generally  take 
the  form  of  either  blimting  the  radome  tip,  which  increases  the  tip  area, 
or  placing  a  protective  metal  cap  at  the  tip  of  a  pointed  system.  This 
section  is  included  to  summarize  the  electrical  effects  due  to  the  blunt¬ 
ing  of  the  radome  tip  which  have  been  indicated  in  previous  sections  and 
to  illustrate  the  far-zone  pov/er  pattern  plot  which  the  computer  program 
generates.  Calculations  of  boresight  error,  transmission  efficiency, 
and  euitenna  pattern  distortion  are  presented  for  an  identical  antenna 
system  operating  with  a  pointed  radome  and  with  a  blunted  version  oi 
the  same  radome.  All  calculations  are  for  perpendicular  polarization 
at  the  center  design  frequency.  The  two  radomes  used  in  this  example 
are  the  pyroceram  ogive  radomes  described  earlier  with  and  without 
the  hemispheric  nose  cap. 

Figure  52  shows  the  effects  of  blunting  the  radome  tip  on  the  bore- 
sight  error  characteristics  of  the  antenna- radome  system.  Figuxe  53 
shows  the  relative  on-axis  transmitted  power  calculated  from  the  sum 
pattern  of  the  monopulse  antenna  for  the  two  systems.  The  "no-radome" 
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Fig.  52.  Calculated  effects  on  radome  boresight  error  due 
to  the  blunting  of  the  tip  of  a  pyroceram  ogive 
radome.  Perpendicular  polarization  at  the  design 
frequency. 

case  is  represented  by  the  100%  relative  power  level  at  all  look  angles. 
From  these  figures  it  can  be  seen  i-hat  the  ray  tracing  theory  predicts 
the  pointed-nose  case  to  give  considerably  better  boresight  error  per¬ 
formance  as  well  as  having  about  5%  greater  on-axis  power  transmittion- 

Figures  54,  55,  and  56  show  the  normalized  far-zone  power  patterns 
for  the  antenna  without  ladome,  with  pointed  radome,  and  v/ith  the  blunted 
radome  respectively  for  look  angle  0®  . 
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Fig.  53.  Calculated  transmission  efficiency  for  the  case 
of  the  radome  of  Fig;  52. 
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Fig.  54.  Computer  generated  far-2one  power  pattern  for 
the  antenna  without  radome. 
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Fig.  55.  Computer  generated  far-zone  power  pattern  for  the 
euitenna  in  the  presence  of  a  pointed-nose  ogive 
radome.  Perpendicular  polarization  at  the  design 
frequency . 
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Fig.  56,  Computer  generated  iar-v.one  power  pattern  for  the 
blunled-nose  version  of  the  radome  of  Fig.  55. 
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V. 


CONCLUSIONS 


A  completely  computerized  two-dimensional  ray  tracing  analyses 
of  radome  boresight  era  or  and  antenna  pattern  distortion  has  been 
developed.  Application  of  this  method  to  several  complicated  cintenna- 
radome  problems  has  been  demonstrated  which  shows  the  usefulness  of 
the  method  both  for  the  design  and  the  analysis  of  antenna- radome 
systems.  Several  example  cases  were  calculated  and  compared  with 
experimental  data  and  with  other  calculations.  Agreement  between 
measurements  and  calculations  was  in  general  reasonably  good.  The 
method  was  modified  to  include  the  analysis  of  aperture  blocking  effects 
in  order  to  form  a  basis  for  the  calculation  of  radome  systems  involving 
matcillic  nosecaps.  Calcttlations  and  meas  urements  of  this  case  will  be 
included  in  a  future  report. 

The  computer  progreim  written  for  the  calculation  of  this  two- 
dimensional  method  is  relatively  long  and  involved;  however,  it  has  been 
written  in  such  a  way  as  to  make  its  use  by  others  relatively  simple  and 
convenient. 

The  method  can  be  applied  to  a  wide  range  of  antenna- radome 
problems.  Calculation  of  results  is  extreme' ’y  fast,  thus  making  the 
method  an  economical  approach  to  radome  de  jign  and  modification. 
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APPENDIX  A 


This  appendix  contains  calculated  boresight  error  (BSE)  and 
relative  power  transmission  (RPWR)  for  several  equally-spaced  ray 
selections  at  te?!  look  angles  to  further  demonstrate  the  convergence 
of  the  ray  tracing  solution  which  was  explained  in  Chapter  IV- A.  Tables 
VI  and  VII  show  the  BSE  and  RPWR  as  a  function  of  the  number  of  equally- 
spaced  rays  (NRE)  used  for  the  pyroceram  ogive  radome  with  and  without 
the  hemispheric  nosecap.  Perpendicular  polarization  was  used  through¬ 
out.  BSE  is  in  milliradians  and  RPWR  has  a  value  of  1.0  for  the  no- radome 
case.  Look  angle  (<t>L)  Is  degrees. 
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Pyroceram  Ogive  With  Hemispheric  Nosecap 
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TABLE  VII 

Pyroceram  Ogive  Without  Nosecap 
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APPENDIX  B 


This  appendix  contains  calculated  boresight  error  and  on-axis 
power  triinsmission  for  several  choicef^  of  phase  allowance  (PHD)  and 
amplitude  allowance  (TRD)  parameters  in  the  numerical  integration 
technique  explained  in  Chapter  IV-B.  The  radomes  used  are  the  same 
as  those  of  Appendix  A.  Lock  angle  (<!>l)  in  degrees,  boresight  error 
(BSE)  is  in  milliradians,  PHD  is  in  degrees,  TRD  and  relative  on-axis 
power  (RPWR)  are  normalized  fractions.  The  number  of  subapertures 
(N)  used  refers  to  the  number  of  subsections  of  the  500  point  recon¬ 
structed  aperture  remaining  cifter  averaging  the  aperture  distribution 
functions.  Tables  VIII,  IX,  X,  and  XT  show  the  calculated  BSE,  R’^WR, 
and  N  for  12  choices  of  PHD  and  TRD  for  the  blunted  (with  hemispheric 
nosecap)  and  pointed  (without  nosecap)  rademes  respectively. 
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TABLE  IX 

Ogive  Radome  With  Hemispheric  Nosecap 


I 
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0.20  -0.03  0.994  2  0.503  0.994 


^  TABLE  XI 

Ogive  Radome  Without  Nosecap 


f 
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APPENDIX  C 

INSTRUCTIONS  FOR  USE  OF 
COMPUTER  PROGRAM 


A.  Required  Input  Data 

The  following  table  lists  all  the  required  inputs  for  a  calculation 
of  radome  boresight  error  and/or  antenna  pattern  properties .  All 
data  are  read  in  from  punched  cards.  The  order  of  the  input  of  the 
data  depends  upon  the  calculation  and  therefore  can  be  varied. 
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TABLE  XII 


Variables 

MC 

LL(I) 

PL(I) 

NRE{I) 

POLIZ 

TITLE 

SOURCE 

FREQ 

NUM 

SPAN 

NOS 

N{I) 

SHAPE!  I) 
DCE(I.  ID 


Description  of  Variable _ 

Niunber  of  cases  to  be  calculated 

Number  of  look  angles  for  which  cal¬ 
culations  are  made  fcr  the  I-th  case 
on  MC 

The  I-th  look  angle  on  LL 

The  number  of  equally  spaced  rays 
used  for  the  I-th  case  on  MC 

Polarization:  Either  perpendicular 
or  parallel 

Description  of  Calculation 

Specifies  source  type,  taper,  pre¬ 
sence  of  obstacles,  presence  of 
radome,  source  taper,  etc.  See 
comment  cards  at  beginning  of  pro¬ 
gram  for  complete  description 

Frequency  in  megacycles 

One  less  than  the  number  of  points 
calculated  in  a  partial  pattern 
calculation 

One-half  the  angular  range  over  which 
pattern  points  are  calculated 

Number  of  geometry  sections  used  to 
define  the  radome  shape 

The  number  of  layers  in  the  I-th 
geometry  section  of  NOS 

The  shape  of  the  I-th  section  on  NOS 

The  dielectric  constant  of  the  I-th 
layer  of  the  IL-th  geometry  section 


__  Units _ 

Integer  Number 

Integer  Number 
Decimal  Degrees 

Integer  Number 

Alphabetic 

Alphanumeric 

Alphanumeric 

Decimal 

Integer  Number 

Decimal  Degrees 

Integer  Nixmber 

Integer  Number 
Alphabetic 

Decimal  Number 
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T.ABLE  XII  (Cor.t.) 


Variable 
D(I,  ID 

TD(I,  ID 

XOO,  YOO 

ROO 

XLIN,  YLIN 

PHLIN 

DFATSP 

A 

TAPER 

PHD 

TRD 

NBLOK 

XGB(I) 

STAP(I) 


Description  of  Variable 

The  thickness  of  the  I-th  layer  of  the 
IL-th  geometry  section 

The  loss  tangent  of  the  I-th  layer  of 
the  IL-th  geometry  sec'.ion 

The  coordinates  of  the  center  of  an 
ogive  radome  section 

The  radius  of  an  ogive  section 

Coordinates  of  a  point  on  a  conical 
section  of  a  radome 

Included  angle  of  a  conical  section 
of  a  radome 

The  distance  from  the  coordinate 
axes  to  the  source  aperture  plane 

The  length  of  the  aperture  plane 

Description  of  aperture  amplitude 
taper 

Phase  allowance  used  in  averaging 
the  reconstructed  aperture 

Transmission  allowance  used  in 
averaging  the  reconstructed  aperture 

The  number  of  a  geometry  section 
which  is  an  aperture  block  (metal). 

If  none  NBLOK  =  0 

The  X- coordinate  of  a  boundary  be¬ 
tween  two  geometry  sections,  I  and 
I+l 

The  amplitude  of  a  step  used  in  a 
step  +  function  amplitude  taper.  I 
on  MC 


Units 

Decimal  Inches 

Decimal  Number 

Decimal.  Inches 
Decimal  Inches 

Decimal  Inches 

Decimal  Degrees 

Decimal  Inches 
Decimal  Inches 

Alphanumeric 

Decimal  Degrees 

Decimal  Number 

Integer  Number 

Decimal  Inches 

■Decimal  Number 


77 


B .  De s cription  Of  A  Typical  Calculation 


The  SOURCE  card  plays  a  key  role  in  determining  the  calculation 
procedure.  SOUF.CE  contains  twelve  alphanumeric  words  which  control 
various  phases  o',  the  calculation.  Twelve  comparison  words  are  read 
into  the  program  as  data  statements.  By  comparing  the  contents  of  the 
SOURCE  card  tc  the  data  statements  the  desired  subroutines  are  called 
or  the  desired  calculations  aie  made.  Each  word  of  SOURCE  contains 
6  alpheinumeric  characters  or  blanks  (which  are  designated  by  . 

Each  word  is  explained  below . 


Source  (1)  indicates  the  calculation  of  the  monopulse  difference 
pattern  by  FlvlONOP.  The  sum  pattern  is  calculated  otherwise. 

Source  (2)  is  not  used. 


Source  (3)  indicates  the  calculation  of  the  no-radome  case  by 

b  . 


Source  (4)  indicates  the  use  of  an  aperture  pedestal  +  function 
taper  by  the  word  STEP'^’^. 

Source  (5)  indicates  the  absence  of  an  aperture  amplitude  taper 
by  FNOAAT. 

Source  (6)  indicates  the  absence  of  an  aperture  phase  taper 
by  FNPHST. 

Source  (7)  causes  the  program  to  calculate  sidelobe  level  and 
half-power  beamwidth  when  the  word  OPTION  occurs  there. 

Source  (8)  indicates  the  use  of  a  circular  aperture  by  CIRCLE. 

Source  (9)  indicates  the  presence  of  an  aperture  block  by 

'’block. 

Source  (10)  calls  for  a  graph  of  the  far- zone  pattern  by  PLOT^*^. 

Source  (11)  indicates  that  special  x-geometry  boundaries  are  to 
be  read  in  by  the  word  SXGB^^. 

Source  (12)  causes  awrite-out  of  ray-tracing  and  multilayer  cal¬ 
culations  by  WRITER. 
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Similarly  TITLE(l)  indicates  the  calculation  of  bo'esifht  error  by 
the  word  FNULL^.  Otherwise  a  partial  pattern  calculation  is  carried 
out.  All  other  alphanumeric  input  cards  are  similarly  used,  e.g.,  the 
words  ^OGIVE  and  ^CONZ^  on  the  shape  card  cause  the  program  to 
call  either  the  ogive  or  cone  subroutine  respectively,  to  calculate  the 

geometric  constants  associated  with  a  radome  section.  One  alpha¬ 
numeric  word  enables  the  taper  subroutine  to  utilize  c,ay  one  of  a  set  of 
pre-programmed  aperture  tapers.  The  TAPER  card  specifies  this 
word. 


Thus  for  a  typical  calculation  the  number  of  different  cases  to  be 
run  and  the  various  differences  between  cases  are  determined.  Data 
which  does  not  chsinge  from  case  to  case  is  read  in  between  the  following 
two  specific  cards: 

If(L.GT.l)  Go  to  937 

Read . 

937  Continue 

This  prevents  unnecessary  duplication  of  computer  reading  time. 

Two  or  more  completely  different  sets  of  calculations,  such  as 
the  boresight  error  curves  fcr  differing  radomes,  may  be  carried  out 
on  one  run  by  utilizing  the  cards 

DO  949  lY  =  1,4 


Complete  program  ! 


949  Continue 

as  shown.  These  cards  are  in  the  program  permanently  so  that  runs 
may  be  combined  without  re- compiling  the  program  (which  consumes 
about  1/2  minute).  The  number  4  is  completely  arbitrary.  These  cards 
cause  an  abnormal  termination  of  the  c  mputer  and  consequent  error 
message  vvhen  less  than  4  runs  are  made,  however  this  is  no  problem 
because  the  calculations  are  finished  when  the  termination  occurs. 

C.  Sample  Data  Lists 

Two  example  calculations,  one  of  boresight  error  and  one  of  antenna 
pattern  parameters,  are  discussed  to  illustrate  the  use  of  the  program. 
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1 .  Boresight  error 


This  example  uses  an  ogive  radome  having  a  hemispheric  nosecap. 
Since  the  calculation  requires  separate  equations  for  the  geometry  of 
the  ogive  walls  on  the  upper  and  lower  sections,  the  cap  is  also  divided 
into  two  sections  to  maintain  symmetry  in  the  geometrical  description 
of  the  radome.  This  results  in  a  four-section  radome  as  shown  in  the 
sketch  below. 


Wall  construction  in  all  four  sections  is  a  constant  A- sandwich.  Five 
cases  are  to  be  calculated  at  four  look  angles.  Refering  to  Table  XIII 
case  MC  =  1  is  the  "no-radome"  case  which  is  calculated  for  reference. 
Since  the  pattern  without  radome  is  independent  of  look  angle,  PL  =  0 
is  used.  Cases  MC  =  2,3  correspond  to  perpendicular,  parallel  polari¬ 
zation  calculations  for  the  same  radome.  Cases  4,5  correspond  to 
perpendicular,  parallel  calculations  for  a  similar  radome  but  with  a 
different  core  dielectric  constant.  All  data  of  Table  XIII  are  explained 
on  the  right  hand  margin  of  the  table.  Reference  to  Table  XII  will  define 
the  variables . 

2 .  Antenna  pattern  parameters 


This  example  uses  an  ogive  radome  having  a  hemispheric  nosecap. 
The  aft  portion  of  the  radome  is  conical  for  fairi:ig  purposes.  Thus 
there  are  6  geometry  sections,  3  on  each  side  of  the  coordinate  aj:es. 
Eight  cases  are  calculated  corresponding  to  four  different  aperture  dis¬ 
tributions  with  and  without  radome.  The  far- zone  pattern  is  calculated 
and  plotted  for  the  angular  range  of  +  45*  (span)  about  the  beam  axis. 
Sidelobe  level  and  half-power  beamwidih  with  and  without  radome  are  the 
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quantities  being  calculated  as  a  function  of  source  distribution.  As  can 
be  seen  from  Table  XIV  the  quantity  of  data  read  in  after  the  initial  run 
(MC  =  1)  is  minimal.  This  is  typical. 
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TABLE  XIII 


VARIABLE 


5 

1 

10.0 

100 

3, 

1. 


.1 


MC 

LL 

NBLOK 

A 

NUM 

PHD.TRD 

SPAN 


501  501 

501  501  501 

NRE 

FNULL  CALCULATION 

TITLE 

9300. 

FREQ 

4 

NOS 

3  3 

3  3 

N 

OGIVE  OGIVE  OGIVE  OGIVE 

SHAPE 

FMONOPULSE 

NO  RADOME  FNCAATFNPHST 

SXGB 

SOURCE 

-2.  DO 

-  2.0D0  3.0D0 

2.0D0 

-2.  DO 

SXGB 

..  -2. DO 

69.0DC  80.0D0 

XOO,YOO,ROO 

24.  DC 

G.DU  2.0D0 

xoo.Yoo.noo 

.  24. DO 

O.DO  2.0D0 

XOO.YOO.ROO 

-2.  DO 

69.0D0  80.0DO 

XCO,YOO,ROO 

PERPENDICULAR 

POLARIZATION 

0. 

t 

1 

PL 

8.9 

2.8  8.9 

DCE 

.040 

.192  .040 

D 

.0003 

.0057  .0003 

TD 

8.9 

2.8  8.9 

DCE 

.040 

.192  .040 

D 

.0003 

.0057  .0003 

TD 

8.9 

2.8  8.9 

DCE 

.040 

.192  .U40 

D 

.0003 

.0057  .0003 

TD 

8.9 

2.8  8.9 

DCE 

.040 

.192  .040 

D 

.0003 

.0057  .0003 

J 

f 

TD 

4 

1 

NOS 

3  3 

3  3 

N 

OGIVE  OGIVE  OGIVE  OGIVE 

SHAPE 

FMONOPULSE 

FNOAATFNPHST 

SXGB 

SOURCE 

-2. DO 

2.0D0  3.0D0 

2.0D0 

-2,  DO 

SXGB 

-2,  DO 

69.0D0  8O.OD0 

XOO.YOO.ROO 

24.  DO 

O.nO  2,0D0 

MC 

=2 

XOO,YOO,ROO 

24.  DO 

O.DO  2.0D0 

XOO.YOO.ROO 

-2.  DO 

69. ODD  8O.OD0 

XOO.YOO.ROO 

PERPENDICULAR 

POLARIZATION 

0. 

8.  16.  24. 

PL 

8.9 

2.8  8.9 

DCE 

.040 

.192  .040 

D 

.0003 

.0057  .0003 

TD 
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r 

s 


[ 


8.9 

2.8 

8.9 

.040 

.192 

.040 

.0003 

.0057 

.0003 

8.9 

2.8 

8.9 

.040 

.192 

.040 

.0003 

.0057 

.0003 

8.9 

2.8 

8.9 

.040 

.192 

.040 

.0003 

.0057 

.0003 

PARALLEL 

0. 

8.  16. 

24. 

PENPENDICULAR 

0. 

8.  16. 

24. 

8.9 

2.4 

8.9 

.040 

.192 

.040 

.0003 

.0053 

.0003 

8.9 

2.4 

8.9 

.040 

.192 

.040 

.0003 

.0053 

.0003 

8.9 

2.4 

8.9 

.040 

.192 

.040 

.0003 

.0053 

.0003 

8.9 

2.4 

8.9 

.040 

.192 

.040 

.0003 

.0053 

.0003 

PARALLEL 

0. 

• 

CO 

24. 

TABLE  XIII  (Cont.) 

VARIABLE 

DCE 

D 

TD 

DCE 

D 

TD 

DCE 

D 

Y  TD 

POLARIZATION 

MC=3  PL _ 

;  I  POLARIZATION 
PL 
DCE 
D 


MC=4  TD 


HC=5  PL 
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TABI.F  XIV 


VARIABLE 


8 

11111111 

2. 7500 
10.0 
100 

3.  ,1 

45. 

C0S6 

501  501  501  501  501  501  501  501 


FPARTIAL  PATTERN  DEG 

6 

111111 
CONE  OGIVE  OGIVE  OGIVE  OGIVE  CONE 
O.DO  7.0D0  -.500 

.00  -69.000  80.000 

6.5000  O.DO  2.000 

6.5000  O.DO  2.000 

.00  69.000  80.000 

O.DO  7.000  .500 

PERPENDICULAR 

0. 

5.5 

,297 


5.5 

.297 

5.5 

.297 

5.5 

.297 

5.5 

.297 

5.5 

.:^97 

9600. 


FLINE  SOURCE  NO  RASTER 

FNFHS'^OPTION 

PLOT  SXGB 

-3.00 

.00 

2.760D1 

3.000001 

2.76001 

.00 

-3.  DO 

0,  0.  .25 

.025 

,05  .05  .071 

.071 

aiNE  SOURCE 

FNPKSTOPION 

PLOT  SXGB 

-3.00 

.00 

2.76001 

3.0000D1 

2.76001 

.00 

-3.00 

aiNE  SOURCE  NO  RA 

FNP-i'SrOPTION 

PLOT 

FLINE  SOURCE 

FNPHSTOPTICN 

PLOT 

FLINE  SOURCE  NO  RA 

FNPHST0P7I0N 

PLOT 

FLINE  SOURCE 

FNPHSTCPTION 

PLOT 

FLINE  SOURCE  NO  RA 

FNPHSTOPTION 

PLOT 

FLINE  SOURCE 

FNPHSTOPTION 

PLOT 

MC 

LL 

NBLOK 

DFATSP 

A 

NUM 

r;iD,TRO 

SPAN 

TAPER 

NRE 

TITLE 

NOS 

N 

SHAPE 

XLIN.YLIN.PHLIN 

XOO,YOO,ROO 

XOO.YOO.ROO 

X00,Y00,R00 

X00,Y00,R00 

XLIN.YLIN.PHLIN 

POLARIZATION 

PL 

DCE 

D 

TO 

DCE 

0 

.TD 

DCE 

0 


MC=1  TO 


f  SXGB 
MC=3  SOURCE 
MC=4  SOURCE 
MC=5  SOURCE 
MC=6  SOURCE 
MC=7  SOURCE 
MC=8  source; 
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D.  Computer  Programs 


HA  Pi 


DATI. 


h9l)') 


C  AIL  CALCUL  AT  "vinw'  Ai.SDM  1 M.;  i^AUDK!-  1  L,  AYi'-N't.  T  '  I C  Ai'l'Jl  r,(-  AX! 

C  AA,  Dfi,cc,,  i)i  ,  I  r- ,  r  f  A'M'  Kill  Detl'pwr;  C')\’ST\:JTS  I-Iv^  a  aiviv  ‘A):'  l;. 

C  AF  N  ll:\'GrH  I'F  '!Th  SUliAPr  H  FURL . 

C  CF  M  PiFA'-.i  CORRiCUON  OUt  TO  THl  ‘mi  SURAPCRTUF’f.  LOCATiC!.\  A'../.Y 

C  FROM  IMF  ORIG'IM. 

C  P  IS  THF  P/.rrFRM  A’iGLE  PHI.  F’FASURF.D  Pi  THT  X-Y  PLANE 

C  PL  1  IS  niL  1  ni  LHO:'  ANGLE  I  DiIGREFS. 

C  AF  IS  Tilt:  FAil  FIFI.l)  AMPI.ITUDF:. 

C  0  IS  IHE  ruT^^ur  PARA'-'FTFR  for  THi-  graph  SUBROllTIOF. 

C  RPL  I  IS  r-HI  LOM-:.  ANGLE  IN  ‘'ADIANS. 

C  FIPI)  IS  TH!  I.NSFRTIO'i  PHASF  OFLAY  FOR  A  RAY, 

C  PH  IS  THb  iniAL  PHASE  TFR*^*  FOR  A  (.  IVEN  S')R  APFR  TURF . 

C  TI  IS  THF  angle  or  INCIDtNCE  FOR  A  RAY. 

c  aat  apfrtiAvF  amplitude  taper. 

C  TITLE  INGICAIFS  WHFTUFR  10  CALCULATE  THP  NULL  UR  A  PARTI  A!  PAIIF  . 

C  BIS  THE  ANGLE  USED  IN  GRAPH  TO  PLOT  THE  PHI  AXIS, 

C  INDICAIFS  THF  DIVISION  HEI^EF,-  THE  TWO  HALVES  OF  THF  TURi'. , 

C  MH  IS  THE  MUNBER  OF  THAT  SUBAPFRTURE  I  D I A  T  F  L  Y  PFLOW  TH"  X-AXIS. 

C  YC  IS  THE  Y-CUORD  OF  THF  CENTER  OF  A  SUCAPFIRTURE  AT  LULiK  A,.r,i,--  o. 

C  PHST  IS  THF  PHASE  TAPER. 

C  DCF  IS  TliE  REAL  PART  OF  THE  RELATIVE  PERMITIVITY. 

C  OCE  IfL  IS  THE  DIELECTRIC  CONSTANT  FOR  THE  ITH  LAYEk  OF  THE  [-TH 

C  CEORETRY  SECTION. 

C  TO  1,1.  IS  THT  LOSS  I  ANGF  'IT  FuR  ITH  LAY'IR.LTH  Gi-ORETRY  S''CTI'‘;. 

C  0  I.L  IS  THl  THICK.NESS  IN  UCHFS  FOR  THE  Ph  LAYER, LTH  ..FGHE  T  <V 

C  TC  IS  THE  AMPLITUDE  TR  ANSY  I  SS  1  ON  COE  F  F  I  C  1  E'''T  FOR  A  RAY. 

C  TR  IS  the  H)TAL  transmission  FACTOR  FOR  A  SUb  APFP.TU'^L . 

C  POLIZ  INDICATES  PERP  OR  PAR  PIU.  A.-l  1  Z  A  T  ION  FOR  VULTILAYFi’.  CALC  . 

C  XGB  IS  THE  X-COORDINATE  OF  A  RADOi-'t  GFOoriRY  D,I  SCMr.'T  1  vij  1  I Y  . 

C  YGP  IS  THE  Y-CnOROlNATF  OF  A  RADOHF  GEO'-' r  TRY  D  I  SFO'NT  P.'U  I  i  Y  . 

C  NRE  IS  THF  NUMMFR  QF  FOUALLY  SPACED  '‘AYS  USED 

C  NKE  IS  CHOSEN  UNEVEN  TO  GET  A  RAY  THROUGH  THE  ORIGIN, 

C  T  IS  THF  ANGLF  OF  P'CICE^!CE  FROM  THF  CENTER  fp  A  SUBAPERTja-  TO 

C  THE  RAOOME  INNER  WALL. 

C  LL  INDICATES  THF  TOTAL  NUMBER  OF  DIFFERENT  LOOK  A’^GLi'S  EXA‘;MFD. 

C  XOO,YOO,ROO  INDICA1E  THE  CENT  ER ,  RAD  I  US  FOR  AM  OGIVi-  SEC11(l^;  M 

C  THE  X,Y  PLANE. 

C  YO  IS  THE  MAXIMUM  MAGNITUDE  OF  Y-COO-’’ f)  I N \T  L  FOR  A  SUrAPF.RIURt  AT 

C  LOOK  ANGLE  0.0  DEGREtS. 

C  MC  IS  THF  TOTAL  NUr'.BER  OF  CASLS  TO  0[  RUN 

C  A  TOTAL  SOURCE  APERTURE  LENGTH. 

C  NOS  IS  THE  NUF-BER  OF  RAOCMF  STCTIONS  HAVING  Oir-FF'^FNT  oEOfLTRY. 

C  FRFO  IS  THF  F»COUr-NGY  IN  MEGACYCLLS. 

C  ANFA  IS  THP  number  OF  EOUAL-LENGTH  SUBAP^PTURES. 

C  NUM  IS  ONI.  LESS  THAN  THE  MIMLER  OF  FIELD  POINTS  TO  BF  CALCULATED 

C  IN  A  PARTIAL  PATTcRN  TYPE  CALCULAIIQN. 

C  FJ  IS  THF  COMPPX  NUMBER  J 

C  SPAN  IS  THF  ANGULAR  RANGE  DF  CALCULATION  ABOUT  THF  PAH  AXIS. 
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ooor;  oooooooooooooooooooooocTOooooooor* 


MAIN 


OATP 


691^)0 


C  I’SKM  IS  llif  I’CT'-FSIGHT  CR^T’O  IN  M I L  L  I  ^  A()  I  ANS 

IS  IS  nil  SLufM-  or  nip;  TANP;r..i  ru  mK  HAufiMt  cu'^vf  at  p*  x,y  . 

TL  IS  TPIti  UF  TUP-  RAY  AT  TIIF  SA’^^  POINT  P  X,Y 

SI.L  IS  SinP-l/'(’,E  LEVEL  IN  PP=R  CLMT. 

SLOE  IS  SKMiU'I.F  LFVFL  IN  OR 
Ml  IS  THE  ’.'jy.BFK  OK  SUBAP’FRTURFS 

I  scale  is  a  SPRFAnr.T,  FACT'TR  FDR  GRAPH  WHEN  FEU  POINTS  APF  USED. 
RPWR  IS  TUP-  POWFR  AT  LCAM  MAXIMUM  RFLATIVPi  TO  THL  VO  RADOf’t-  CASE. 

Li.x  IS  THL  p;H-si',n  VALUE  nr  the  look  angle  index  i  in  main. 

THE  TLIM.  VARIABLE  PIU-VCNIS  COMPUTER  ROUNDOFF  rRRfJP.. 

FNK[>  IS  IHC  I’DUF-'  JORMAL  1 7  INC.  TERM. 

TAPER  l.NDICAIESI  IN  A  ONIf  WORD  HOLLERITH  FIELD)  THF  SOURCE  TAPER. 
OP’.ST  IS  1/P  THF  LENGIII  OF  A  SY.MMi  TR ICALLY  PLACED  AP6RTIJP.P.  I’LOCP;. 
SOUPCL-Il)  IVniCATrS  the  calculation  of  a  MONOPULSF  DI  Ff-Ef'i:  \ci- 
PATTFRN  DY  ’MONO*  THE  SUM  PATTEPfl  IS  CALCULATED  nTHrKWlS‘^. 
SnURCE(3)  INDICATES  CALCUI.AIION  OF  THE  Nf)  RADOME  CASE  hY  '  :\U  • 
SOURCtlA)  INDICATES  TUP-  USE  UF  AN  APERTURE  PLDESTAL  +  FUNCT  I U-N  TAPER 
PlY  THE  Kt;.(CO  'STEP*. 

SOURCE  (S)  INOICATPS  THF  ABSENCE  UF  AN  AM'^LITUDE  TAPER  UY  ’ENEjAAT'. 

SOiiRCF(A)  luniCATLS  THE  APiSE.’sCh  UF  AN  PHASE  TAPER  BY  'FNPHST*. 

SOURCE  7  CALI.S  GP' APll,  SLL ,  HPBW  WHEN  THE  WORD  OPTION  OCCURS  THtRF. 
SOURCr.  (fl)  INDICATES  THE  USE  i<F  A  CIRCULAR  APFRTUNfr  BY  'CI  ’CLE*. 

Sf)URCi:(9)  INDICATFS  THE  PR(  SI  NCE  OF  AN  APF.UURE  BLOCK  BY  'PanCK', 

SHURCI:  (  ID)  C  ALLS  FOR  A  GRAPH  Of  IP'P-  FA!<-riELr;  bY  ‘PLOT*. 

SOURCF(ll)  r.'OlCATES  THAT  SPECIFIC  X-(;i  OMfi  TRY  HOU‘'DARIFS  ARE  TO  c.E 
REAP)  IN  f.Y  'SXGB'.  OTHERWISE  PROGRAf^  CALCULATES  THi-  XGP^'S. 
S0URCC112)  CAUSES  A  WRITF-OUT  OF  ,<AY-TP.ACF  A;^^n  MULTILAYER.  CALCUL-U 
IONS  IF  THE  UORD  WRITER  OCCU.RS  TI'CRF. 

THE  DO  9-^9  LOOP  ALLOWS  THE  ENTIRE  PROGRAM  TO  PC  ''L-CYCLLD. 

RAF  ,  RCF  ,  .<rRi  Kf’H  arc  TEMPORARY  I.OCATIONS  FOR  P  FCONS  T  RUC  fFD  APr.P.TUr’.F 
C  0  L-  E  F  1  C  I  N  T  S  A  f  = ,  C  F  ,  T  R  ,  P  0  . 

DFATST  =  DISTANCE  FROM  AXIS  TO  SOURCE  PLANF. 

THE  SIOELCI'E  LEVEL  SUBROUriNL  CALCULATES  TI'F  POWf.P'  LP'.VlL  AT  THE  1ST 
LOP<F  LPJCATP^N  OF  THE  fv!G  MADOME  CASE  EVIRY  Tl.ME  IN  ADDITION  TO  T!-- 
MAXIMUM  SIOt-LOf’e  LEVEL.  WHEN  NO  SIDLLOBF  15  DITLCTED  THIS  IS  TAP’-.L'. 
SIDELORT;  LEVEL. 

DFATST  =  DlSTAPfCF  FROM  AXIS  TO  SOURCE  PLANE. 


1)01  IKLF  PR.f-C  1  S  I  ON  AA.BR  ,CC  f  00,1.  e  fEE,  XGP5,  YGI^.P  I  ,  BBH  ,r.CC,  DI  SC  ,  DYDX,  X, 
lY,  TS,TI  ,CL  ,CLPf  SL,  SL?,C?l.,S?Lf  Tl  ,  iU>L  ,  D  1  SCU  t  XOC ,  YUD ,  Rf.O 
1  ,  XA,YA,DEATSP,XLIN,  Yl.  IN.PHLliN 

COMMON  T  (  GO  I  ,  I  ?  )  ,  SOURCt.  (  I  2  )  t  P  I  ,  RADE  5,01;  GR  AD,  M,1  ,  NOS  ,  I  X  ,  M.C  ,  A  ,  YC  (  SO  1  ) 
I  ,  XA  (  SOI  )  ,  Y/»(  son  ,Dr  ATSP,  XLI  N,  YLIO,  PHL  IN, 

1  P'llSTCSOl  )  ,AAT(SOl)  ,  T  A  PEP’. ,  fPBLOK  I  ,  NHL  0K2  ,  YO  (  SO  1  )  ,  VCX 

CnP/MO.N  /TPR/  Y,STAP(2S) 

COMMON  /RR!/  ,RPL  (  2  S  )  ,  XGR  (  1  ?  )  ,  T L  ,  SL  ,  S  ?L  ,  CL  ?  ,  S L? ,  XOO,  YOO ,  ROf.  ,CL  , 

1  AA(  12  )  ,HP.(  12)  ,CC(  12)  ,DD(  12  )  ,  EF  {  12  )  ,FF(  12) 
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MA  I  .‘'j 


A  T  l_ 


f m  ^  0 


i.an'ca,  at { ‘j(;i  )  ,r,i 
1\''?  ,  L  I  X,. 

COKM()\  /•  JKV/  \'(  12  )  ,  i)(  12,  I?)  ,OCL  n?,  I  ?  )  ,Ti;(  1?,  1  ;^ }  ,F  IP::(  bOl  )  ,  TC<  “jO 
l,FRFg,POL17(  n 

niMEgsiri  LL(2S),  IPl'SOn,  7  I  TLCi  12)  ,  .l 

I  1  ),  AH  (  501  ),  Pill  F.oi  ),  AC  (‘30  n  ,P  (  1  2  >  tAL  (  ^on  ,  SHAPE  (  1  2  )  ,  '’HDM  bC  I  ) 
niMFNSIO'i  ".AF  (  501  )  ,  HCF  ( ‘301  )  ,  RPH  (  OO  I  )  ,  PTR  (  50  !) 

1  FORMAT! lOFP.O) 

2  r  URHAT  I  I  5  ) 

3  FORMA  r  (  \  15) 

4  FORMAT  (F  10. 3, hH  SFCOfinS) 

5  FORM, AUF  10,3, HH  yiNUTES) 

9  FORMAT ( 601 3.8 ) 

10  FORMAT! lo I  5) 

11  FORMAl  (8F10,6) 

1?  FORMAT! 301 0.6) 

13  FORMAT  ! IS.FIO.A, 15) 

14  FORMAT!  I  Hi  LOOK  A'lCL  t ,  F 10. 3  ) 

15  FURMAT(3H  MC,15/)liJ  WA  VF  Lt  .'iO  TH ,  F  1  0. 6/ ?9H  I'.’O,  OF  APCRTURF  SUHOlV! 
inNS,I5/28H  NO.  OF  LOOK  ANGLES  L  X  AM  I  NOO  ,  ! ‘j  /  /  ) 

16  Ff!RMAT!22H  LOOK  AR’GLL  IN  RAOl  ANS  ,  r  10 . 6  ) 

17  FORMAT (F 10.2) 

18  FORMAT! OF  ).6) 

19  FORMAT  (H;l) 

23  r-0kMAT!6i;  12.4  ) 

40  FORMAT ( F 1 ?  . 2 , F 1 1 . A , F 1 1 . 4  ) 

42  FORMAT! 13F6,?} 

45  FORM.ATIFI  i.5,F9.5) 

100  FORMAT!/!?! 10.4) ) 

140  FORMAT!  )2H  SPAIi  A’T.l.F  ,F10,4) 

141  FORMAT! 19H  PfOCSTAL  HEIGHT  ,F10.5) 

150  f  ORMAT! l?(  A4,2X)  ) 

151  F0RMAT(3( A4,2X) ,F4. I ,6(A4,2X) r 14,a4,2X) 

152  FORMAT ! 4X, 17 Aa) 

COMPLEX  FJ,F  !‘)C1  )  ,1:3 
DATA  STEP/4HSTFP/ 

OATA  PARALL/4il  PAR  / 

OATA  WRI  TLH/'fHKRI  T/ 

DATA  SXGIJ/4HSXGB/ 

DATA  PLnr/4HPL'^T/ 

DATA  FMULL/AHFnUL  / 

DATA  FMOX.'lP/'tHFMOM  / 

DATA  FNnHM)/4H  NO  / 

DATA  OGIV‘//.M  Of.I  / 

DATA  C0Nt/4H  COM  / 

OATA  nPTI0;-i/4H0PTi  / 

PI =3.  I4l 592053 58P 7932 
OtGPAD=PI/ IHU. 
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MAIM 


DA7f 


C9  1';9 


RADtf'.  -  I  BO.  /(>  1 

(•>)=(0.  ,  1.  ) 

DO  9A9  1Y=^1,2 
Kn2=i 

Rt-ADIS,?)  MC 
RFAD(r>,ii)  PH!), Tan 
PHO=PHO*r;rr,RAn 
RHADIt.,?)  ILL!  1),  1^1  ,MC) 

DO  300 
MGl=-0 
MMMM-0 
RP'..'R=1  . 

LT=LL(L) 

MCX  =  L 

ir (L.GT. I )  GO  ro  937 
RtADIG,.?)  M(3L'TKl,MMnK2 
RF.AD{!S,1?)  OFATSP 
READ! 5. I )  A 
FMRP=A 

READ(!i,?)  MUM. 

NUM2=MUK+ i 
TNUM  -NU.M 
Tf«IUK2  =  Tinj  -VP. 

RCAD{S,17)  SPA^! 

RSPA.\’  =  S'’AM  =  Di..''.RAG 

read(5,i‘jO)  taper 

RCA0(‘>,2)  (.MRE  (  I  )  ,  I  =  1,MC  J 
RtAf)Iii,190)  (TITLE  (I),I  =  l,l2) 

READ  (*5, 2)  MOS 
NNOS^l+.NOS 

READCitlO)  {M(I),I  =  l,MnSI 
REAO(5,1‘jO)  IF.HAPh  (  I  )  ,  1  =  1,  MOS) 
no  121  1=1, NOS 
IX=I 

IF  (SHAPE  (  IX).  r  O'.  OGIVE  )  GO  TO  155 
IFISHAPbi  IX).l.O«CCMS)  GO  TO  157 

RE  AD  (5, 9)  AAdX)  ,rR(  IX)  ,CC(  IX)  .DDI  IX)  ,EEl  IX)  ,1  F{  IX) 
GO  TO  121 

155  PLAD(5,i2)  xon,YQn,Hnn 
CALL  GIVE 
GO  TO  121 

157  P.EAD{5, 12  )  >L  IM,  YLIM,PHLIM 
CALL  CONIC 
121  CONTINUE 

RFAD(5,150)  (POL  I Z{ I ) , 1  =  1 ,3) 

READ (5,1)  ( ’E( i ) , I  =  l,LT) 

IT (POLI? ( 1 )  .rO.PARALL)  GO  TO  909 
DO  120  IL=1,N0S 
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MOO  i 


y.fl  I M 


n/.  i  u 


f-.<>  1  ‘I'l 


N'L=M(  II.  ) 

Rt  AIK  I  1  )  ( '.C-t  (  I  ,  I  U  ,  I  -  I  ,\L  ) 

RFAIK'i,!!)  (  IK  I  ,  :L  )  ,  1  =  1  ,.',1  5 
120  RFAnrKll)  <  Ton  ,  iL  )  ,  I  =l  ,‘!L) 

909  CONTINUF 

WRITh (6, 19) 

RtAO  (3,11)  FRFO 
Y=1 .  18031 AF+A/FRl 0 
W=2.*P1/Y 

937  CONTINUE 

RFA(K‘3,  130)  ISnuRCF  (  J)  ,  J  =  1  ,12) 

IF  (  SOURCI- (  1  1  )  ,  r ’.KSXGfK  i'lAI  (5,9)  (  <05  {  I  )  ,  )  =  1  ,  N’N'OS  ) 

ir  (  SnURCL- (  A  )  ,f  0.  S  TFP)  ni:AD(3,42)  (  S  T  AIM  I  )  ,  1  =  1 ,  NC  ) 

I F  (  SUIIRC  L  (  3  )  .  F  •;> .  F  MGR  AD  )  '*0  1  -  1 
iriTlTLFI 1  ),t0.ruUlL)  GO  TO  938 

WRl  TF  (6, 13  1  )  (T  ITLF  {  I  )  ,  1  =  1,3)  ,SPA\‘,  {  r  ITLF(  I  )  ,  1=3, 10)  ,’;uw?,  1  ill  -  (  1/ 
GO  TO  939 

938  WPITt{6,15?)  TITLF 

939  WRITE(6,15C)  ( SOURCE { 1 J ), IJ=1 , 1? ) 

WRITE  (6, 130)  (SHAPE  {  I  )  ,  1  =  1  ,N(;S) 

WRITL{6,150i('nLlZ 

WRITE(6, 152)  TAPER 

IF(STAP{L  )  .GT,0.  )  WRITE(6,  lAl  )STAP(l.  ) 

DO  no  1  =  1  , LI 
AOEA=NREIL)-l 
M1=ANFA 
NR=NRF(L) 

DO  50  K=1,M1 
AF(K)=) ./ANFA 
T1=K 

50  CF  {K)=Tl/AN£A-(  AMCA+1  .  )  /  (2.*A\'FA) 

LLX=I 

RPL( I )=DFGRAD«PL(  I  ) 

CL  =  DCCS(RPL(  I  )  ) 

I  )  =  ANEA/?. 

SL  =  DSI!'KR‘’L(  I  )  ) 

CL2  =  DC0S(RPL I  1 ) )s»2 
SL2=0SIN(P.Pl.  (  1)  )  **2 
S2L  =  0SIN(?.«Rf’L(  I  )  ) 

TL=SL/CL 
CALL  SCLOKl 
CALL  iMKRRl 
941  T1VF=PCL0KI( 1. ) 

TlPFM=TIMF/60. 

WRITL(6, 1001  ■;  TIME 
WRITE  (0,3)  Tn^EM 

1001  rORMAK IX, 17HT IME  IN  RAY  TR ACP , F 1 0 . 3 , PH  SECONDS) 

IF(SnUP.CE(  3).EP'.i  S'ORAD)  GO  TO  337 
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M  A  I  N 


DATf- 


CALL  iiCl.OKl 
CALL 

r  I:'C----<CLl)K  1(1.) 

T  IMrv=TI,‘',F/60. 

Wl- 1  TCI  6,  1002)  TIMF 

100?  f-n^flAM  IX,  lOtiTIMf;  IN'  ^.U'  T!LAYE:R,F10.'5,8II  SCC0'''DS) 

WRIT!- <6,')) 
bS?  CALL  SCLOAl 
CALL  ATAl>i'K 
TI!'i-  =  RCLn!' n  1.) 

TINry^Tl!‘.>-/60. 

WKI rr ( 6, 100^ )  TIMF 
WRITL-16,b)  Til'FM 

1003  r OR  MAT ( IX, 1  OUT  I  ML  TO  C ALC . T APCH S , T 1 0 . 3 , RH  SECONDS) 

IFISOORCT  (  1)  .NQ.FNn.lAD)  GO  TO  609 
00  SbB  NA=l,Vl 
F'PD(NA)=0.0 
TC  NA)=1. 

558  Cn-niNOF 
609  C(if  'NJF: 

no  '  11=1, XI 

TR(  .  •  )=TC( I  1 )«AAT( 1 1  ) 

PH{  ,  )=Fr'n(  I  1  )<DEGRAD+PHST(  II  ) 

710  CCj'iT  ,UF 

KPITi  (6,410)  P!I(),TRD 

410  FORI-i.M  (2X,?6.(ALLnwF0  PHASE  DIFFtKFNCE  F  J  0.6/?X,  33HALLOWED  TRANS. C 
lOFFF.  niFI  CKt  HCf;  F10.6) 

CALL  SCLUKl 
J1  =  0 
J2  =  0 

405  JOv-KMlI) 

J  =  JO 

1F(J2.F0.0)  G"*  TO  411 
J5=Jl 

jn=jn+i 

j= j + 1 

41?  jn=jri-i 

TR]=T'U  JO) 

Pill  =  PH{  Ju  ) 

413  J=J-l 

1F( J.FQ.l )  GO  TO  414 
DTR  =  TIU  J)-TR1 
nPH  =  Pll(  J  )-('Hl 

1  F  {  A  B  S  {  D  T  R  )  -  I  I ) )  415,415,414 
415  IF(  A3S(DPH)-P1I())  413,414,414 

414  PH?=0. 

TR2^0. 

AF?=0. 
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MAIN 


DATK 


6019-7 


Ch?  =  (). 

no  ^16  ji''=j,jr) 

CF?=CK?+CF( JIM 
AF?=AF?+AF (JIM 
TR?  =  rR2+II^  (  JPl 
A16  PH?  =  r'H?  +  PH(  JP) 

Jl=Jl+l 
T3  =  Ji}-J-H 
RAF ( J1  )=AF? 

RCF( J1  )=CF2/T3 
RPH( J1 )=PH2/T3 
RTR( J1 )=rP2/I3 
jn=j 

[F( J.NE.  I  )  no  TO  A  12 
GO  TO  A18 
All  J?=l 
AOA  J0=JU+1 

TRl=TR( JO) 

PH1=PH{ JO) 

AOl  J=J<-1 

IF( J.FO.Ml  )  GO  TO  A02 
nTR=TR( J)-IRl 
DPH=PH{ J)-PHl 

ir( A8S(nTR  )-T^n)  A07,A07,A02 
407  iF(  A8S(ni>H)-;>HD)  A01«A02,A02 

402  PH2=0. 

TR2=0, 

AF2=0. 

CF2=0. 

00  A03  jp-jn,j 
CF2=CF2+CF ( JP) 

AF2  =  AF2-!AF(  JP) 

TR2  =  TR?  +  T'-’  (  JP) 

403  PH2  =  PH2  +  P1I(  JP) 

Jl=Jl+l 

T3  =  J-JUH 
RAF{ Jl)=Ar? 

RCF( Jl)=CF2/T3 
RPH( Jl )=PH2/T3 
RTRI J1  )  =  TR2/T3 
JO^J 

IF(Ml-J)  405,405,404 
418  Ml=Jl 

Ml.  =  Ml-J5 
MM( I )=ML 
J4=Jl+l 

WPIT£{6,‘i2l  )  JlfMil  ,  J4,ML,J5 

421  F0RMAT(2X,  *  Jl  =  ' 15/'  Ml  •15/»  J4  '15/'  ML  *15/'  J5  M5//) 
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MA[\ 


DAT!- 


h<)\  )0 


C  RAF  ,KCI  TFMPOKARY  LOCAIIHWS  FO?  P.F  CtV.STR'JC  I  FD  Al’F^rU'L 

C  COf  Ff  s  AF  ,CF,  rR,PH. 

WP  I  FF  U.,  A06  )  Jl 

406  FORMAT  {  ?X,  O'jHNURHFP  OF  SUP  APbRTURF  S  USED  IN  RFOOUSTRUC  I  fO  APFRTU^F 
1  15) 

DO  419  J=l ,ML 
AF( Jl^KAI ( J4-J) 

CF( J)=RCI ( J4->  ) 

TR(J)=RTk(  J4-.!) 

419  P0( J)=RPh( J4~J  ) 

DO  4P0  J=l,J5 
AF(  JM',L)=RAr(  J) 

CF(  J  +  KL):^PCF(J) 

TK( J+ML)-RTR( J) 

420  PH( J+ML ) =RPH{ J ) 

DO  417  J=l ,M1 

417  P  H  n  G  {  J  )  =  P 1 1  (  J  )  ♦  P,  A  D  F  9 
WRiT£{6,408) 

408  FURMAT(2X,  35HPFCO,\'STRUCTED  APFRTURF  COC  FF  I C  I  EOT  S/4X  ,  2HAF  ,  8  X ,  2HCF  ,  1 
IXfPH-^R.SX.  2hPH//  ) 

WRITE (6, 409)  (AF i J) ,CF ( J ) , 1 R ( J } , PHDG ( J ) , J= I , M I ) 

409  FORMAT(3F10.t),F10.4  ) 

T1ME  =  RCLG;K1  (  1.  ) 

TIME*''-!  li  ■■/bi  . 

KRITl (6, 1009)  TINE 
WRITE  (6, 5)  TIN'EM 

1009  FORMAT!  IX,  36liTIME  TC  CALC.  RECONSTRUCTED  APERTU»  E  ,  F  1 0. 3  ,  EH  SECONDS 

1  ) 

IF{ I .GT. 1 )  GO  TO  943 
94  3  CCJK'TINUF. 

IFIStJURCL  (  1  I.NF-.FyrjNnP)  GO  TO  25 
DO  27  NM=l,ML 
AF (NN)=-AF (NN) 

27  CONTINUE 

25  WRITE(6,l‘j)  ,MC,Y,mi,lT 
WR1TF(6,3)  MMin 
IF  it;  TLf:  (  1  )  .lO.FNOLL)  GO  TO  60 
GO  TO  34 
60  K=0 

CAIL  SCLUKl 
KR-0 

Pil)=  RSPAN 
PI  2 )=-RSP AO 
500  K=K+1 

F{K)  =  (0. ,0.  ) 

on  900  M=l,Ml 

S  =  WeAttSIGMP(K)  )«AF{M)*.5 

IFIS.LT.O,  lO.A.ND.S.GT.-O.lO)  GO  TO  8(jO 
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DATv: 


SS  =  5  I  N{  )/S 
GO  TO  ,H10 
800  SS=l. 

810  ()  =  W«A*CF(M)*SI:N{f>{K)  )-PH(M) 

co=cns{0) 

5nr.sp,'{0) 

Fl=  TR(M)*A*AF(M)*SS 
E?=  El*CO 
F3=  El*SO*f J 
900  E(K)=F(KI+t.  ?  +  r-3 
890  Al;(K)=CAfiS{L(K)  ) 

IF  (K.EO.l.AMO.KR.rC.O)  GO  TO  t>00 

IF (KR.EO. n )Gn  TO  S02 

KR=KR+l 

IF(  AE(  1  )-AF(  2)  )  ‘301,502,503 

501  P(2)=(P( 1)+P(?) )/?. 

K=1 

GO  TO  500 

503  P(l)=(P(l)+P(?))/?. 

:<=o 

GO  TO  500 

502  BNULL=(P( l)+P{2) )/2. 

TlF‘.E  =  RCLnKl  (  1.  ) 

TI  ‘‘■EK  =  TIO'.  /60* 

WRIT£{6, lOOA)  TINE 
WRITE  (0,5)  TI^’L-M 

1004  FORfMTI  IX,  17HTIME  TO  C  ALC  .  NULL  ,  F  10 . 3 , 8H  SFCD.^HS) 
on  59  NJ=1 ,ML 
59  AF(NJ)=A8S(AF(NJ) ) 

34  CALL  SCLOKl 
N'J3=1 

IF( SOURCE! 7) ,EO.nPT ION)  NU3=NUM2 
DO  92  K=1,KU) 

F{K)  =  (0. ,0,  ) 

Tl=!< 

IF(NU3.EQ. I )  GO  TO  510 

P(  U=  RSPAN-(T1-1  .  )  *KSPAN/TNU)'.2 

GO  TO  511 

510  PI  15=0. 

511  DO  90  M=l,fU 
S=W«A«SIfJ(P(  1  )  )*AF(M)*.5 
IFIScLT.O.  lO.ANG.S.GT.-O.lO)  GO  TO  80 
SS=SIRM  SI  /S 

GO  TO  81 
60  SS=1. 

81  0=W*A*CF| V)«SIN(P( 1 ) )-PH{ V) 

C0=CCS(0) 

S0=SIN(O) 
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ei=^  TR(M)<A*  Af-(r')«5S 
C?--  Ll<0 
L-3=  tl*S'>fJ 
90  P (K ) =F (K ) +F2+H3 

prK)=-(P{  1  )^Pin.(  I  )  )«PA0EG 
09  AE«K)=CABS(E{K)  ) 

!r(N'U3.F0.  l.A'ir.5!in9CC(3>.f:0.F\0'-^An)  FNRr'=\F,K  ) 

1  r  (  K>U3.E0. 1 .  A\‘D.  SnURCf.  (  3  )  .NL  ,F\'OR  AO  )  RnW‘>>=  (  Al  (  K  )  /FNRP  )  ««2 
!  F  ( ;<  .  F  0  .  ( "a  J  3  /  ?  )  .  A  n  .  s  r M  J  R  C  r  (  3  >  .  F  0 .  F  N  0  R  A  D  )  F  \’  R  f •'  =  A  F  {  K  ) 

IF(K.F.tj,  {MIJ3/2  J.AMP.S0UkCF(3)..>lF.FMaRAD)  RPWK  =  (  {  K  )  /  F9RP  )  «  <  ? 

92  CONTINUE 

WR17E{6,t>0A)  RPWP 

50^  FnRMAT{//,2X,26HRELAl IVF  POWER  Al  PNULL  ,F10.6) 

TIMF^RCLOKl  (  I .  ) 

TIMEP  =  TI»-<E/60. 

WRITE (6, 1000)  TIME 

1000  FORMAT! IX, 28HTIMfc  TO  CALC. PARTIAL  PATTER N , F  10. 3 , 0H  SECONDS) 

WRITE (6,0)  TIMFM 
33  CONTINUE 

IF(TITLE( 1  ).uO.FNULL)  GO  TO  101 
y.MAX=  0.0 
nMAX=0.0 
DO  200  J=l,NUM2 
ir(AF(.J).GT.X."AX)  GO  TO  2C0 
GO  TO  200 
200  XMAX=AE(J) 

8MAX--^B{  J) 

200  CONTINUE 
101  COfjTINUE 

IF  {  Sf'.URCE  (  1  )  .NE.FMQN'OP)  GO  TO  611 
IF(TI TLE( 1 l.LC.FNULL)  CO  TO  610 
DO  62  J=i,rjU;-'2 

IF  (  AH  J)  .LT  .  AF  (  J-1  )  .AND.  AE  (  J  )  .LE.AF  i  J  +  1  )  .  AND.  APS  (  P  (  J  ) -PL  (  I  )  )  .LT.2. 
1  )  GO  TO  6  3 
GO  TO  62 

63  XMULL=AE(J) 

BMULL=b( J ) 

62  CONTIiNUE 

BSEM:--(bN'JLL-PL{  I  )  )•  DEGRAD*  1  000  . 

GO  TO  601 

610  CONTINUE 
BStM=»NULL*1000. 

P.NULL  =  BNULL«KADEG 

601  WRITt(6,^0?)  PNULL, nSEM 

602  FORMAT(/FH  DNULL  .F10.6,9H  DEGREES, /7H  BSFM  ,F10,6,IAU  MILLIRAP13 
INS.  ) 

611  IF(T1TLE{  1  ).uO.FNULL)  GO  TO  111 
WRITL(6, 19) 
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M  A  I  M 


nATr  =  i<.9 1  -JO 


ir{scu-"Ci  ( f')  .\i  .HPriON)  co  re.  i\? 

If(Si'AIKCf  (  in).!JI:.rLnT)  go  H)  113 
CALL  SCLt'MCl 
ISCA1.E  =  60/NUM 

CALL  GRAPlil  AL-,XKAX,0.  ISCALLl 

TI  MH  =  PCL()IU  {  1.  ) 

TIMCN‘-Tiyr/6{). 

WRnC(6,lU06)  TlMt 
VjRITe(6,G)  TIMF.N 

1006  FORMAT  (  IX,  13HT  IVF.  INI  GRAPH ,  F  1 0 . 3 , 81 1  SFCONDS) 

113  CALL  SCLOKl 

CALL  HPBIMAE,  XMAX,NUM2,  R) 

TI  ME-=RCLOKl  (  1.  ) 

TlMFM=TIRF/60. 

WRITE (6, 1007)  TIME 
WRirf;.6,0)  TIMtM 

1007  FORMAT!  IX,  12HTIME  IN  HPRl. .  F  10. 3 , 8H  SECONDS) 

CALL  SCLOKl 

961  MMMH=l+MMMK 

CALL  SLL!  AFtXMAX,NUK?,  AB,BSLHAX,AL,B,  T.950,  i004,M.Gl  ) 

GO  TO  951 
964  MKMM=0 

950  IF(S00RCt(3).LQ.FNGRAD)  GO  TO  96? 

WRITE (6, 1551) 

1551  FC>RMAT.:///‘THu  POWER  LFVEL  AT  THE  POSITIUN  OF  THL  1ST  :,irL  LOrE 
1  WITH  Ntl  PADDME  IS  CALCULATED  BELOW'//) 

P{ 1)=RSLKAX 
K=l 
KR=l  I 
P(2)=P(l) 

GO  TO  952 

951  K=0 
KR=0 

P(l)=  BSLMAX-DEGRAD 
P(2)=  BSlMAX+DEGRAD 
1500  K=K+l 

952  E(K)=(0.,0. ) 

DO  1900  H=1,M1 

S  =  W*A«SIN(I’(K)  )»AF(M)*.5 
IF(S.LT.0.10-AND.S.GT.-0.1())  GO  TO  1300 
SS=SIN(S)/S 
GO  TO  1810 
1800  SS=i. 

1810  0  =  W*A«^CF(M)*S1N{P(K)  )-PH(M) 

C0=COS!0) 

SQ=SIN(0) 

El=  TR(M) «A« AF (M) *SS 
E2=  tl*CO 
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1900  r- (  K  )  =1- { f' )  ^  ^?^ ':  ’3 
1090  Al.(K)=CACS(I-(K)  ) 

IF  (K.CQ.  1  .  AMO.KR.F-O.O)  GO  TO  1500 
IF  (KK.F.O.  1  I  UiO  TO  1502 
k;i  =  k-^,+  I 

IF  (  AL-  (  1  /  -A^(2)  )  150  3, 1502, 1501 

1501  P(?)  =  (P(  1  )+fM?)  )/?. 

K=1 

PD  =  P{2  )«kAf.'EG 
U'RITL  f6, 1  3  }  K9  ,PD,K 
GO  TO  i50C 

1503  P( 1 )  =  (P( 1 )+P(2)  )/2. 

K  =  0 

PO^PI  1  )*,<A[)Ff;, 

Wfl.rfc(b,l3)  KR  ,PD,K 
GO  TO  1500 

1502  PSL  ={P(  n  +  P(2)  )/2. 

IF(  (P{  D-Pi?)  ).LT..  00001)  AE(2)^AE(  1) 

SLMAX  ^(AC  (  I  )+AF  (2)  )/(2eX'''AX)*100. 

SLDP=20.« ALOGIOI 100. /SLMAX i 

PSLD=P5L*  <AOFG 

WRIT!;  (6,  1  5  50  )PSLU,  SLOB 

1550  FO'JfV.TI  '  LOCATin-J  OF  MAXIMUM  SlOtLOBc  =  SFIO.A,/'  SlUtLOf’L  LLVtlL 
1=  • ,F10.A, •  DB* 

IFIKMMM.EO. 1  )  GO  TO  961 
96?  CONTIMUE 

TIMr;  =  !<CLOM  I  1.  ) 

TiMeM^TIM;:/60. 

WRITE{6, lOOB)  TIME 
WRITE (6, 5)  TIMEM 

1008  f-ORMATI  IX,  I  lOTIME  IN  SLL,F10.2,8H  SEC.PMDS! 

GO  TO  112 
111  CONTINUE 

WR 1 1  F ( 6 , I N  ) 

112  CONTINUE 

iriSUURCFI  1  ).NE.Ff'ONOP)  GO  TO  IlO 
110  CONTINUE 
300  CONTINUE 
949  CONTINUE 
STOP 
END 
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<>  I  V  f 


I  l<\  I  t  =  <>  V  1  ■/'( 


sui’.Riajf )  r.  I  vf- 

nC'Ul'-t.t  c  I  s !  i;:';  aa,  Hh,cr.  ,nn,  ft ,  h  F  ,  Xf."  ,  Yf,i 

1 Y ,  T S  ,  T  1  ,  C  L  ,  C L?  ,  SI.  ,  SL?  , C ?L  ,  S?L  ,  1  L  ,  f<PL  ,  n  I  sen ,  XOn , 

1  ,  XA,  YA,I)F  ATS!',  XL  I  'J ,  YL  ,  f’HL  I  \ 

CnF'MC\-  T(  001 , 1  P)  ,  SUURCLC  1?)  ,  (M  ,  RAOF'.', ,  DEa^AI),  Y  1  IX, 

I  ,XA(  -501  )  ,  YA(  001)  ,  PFATSP,  XI  I N ,  YL  I  o ,  PMI.  IN  , 

1  PHST  (001)  ,  AAT  (  001  )  ,  T  APCR  ,  \'BL. -K  1  ,Nr.LC>K?  ,  Y0(  001 
COMMON  /.Oil/  KPL(20),Xr.B(l/),Tl.,^L,S'’L,,CL?,SL-^, 

1  A  A  (  I  2  )  ,  BO  (  1 2  )  ,  CC  '  1 2  )  ,  00  (  1  2  )  ,  f  (-  (  I  ?  ) 

I ,  ANF  A  ,  Af  (  00  1  )  ,  CF  {  001  )  ,  CM  (  20  )  ,  L  T  ,  PL.  {  ?0  )  , 

INR,N\0S,LIX,C0? 

I  =  IX 


BO.CCC, 

YUn,RlU) 


)  ,MCX 
XOO,  Y  jfl 
,Ft-  {  12) 
NRi.  (2 


'I  SC', 

MC,  A  ,  Yi; 

,0'T,CL 
0)  , 


AA( I  )  =  i. 
BB(  •  )  =  1 . 
CC( I )=0. 


0D( I )=“2. »XQU 
FEl I )=-?.*Y00 

FF(  I  )=--xnn«'xOG+Yon«YOo-wooeRno 

RETURN 

END 


Fig.  58.  Ogive  geometry  subprogram - 
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cnjic 


Dii]  ( 


691  '9') 


SUI'.i'OUl  irv[  CC  --1C 

OniMUt-  pl^CISH  "I  ,  r  I- ,  X^i’ ,  Y  ,  P  1  ,  I'-O'.,  CCC  . 

lY,rS,  I  l,Cl  ,(:i  P,Sl  .SlP.r.PL.SPL.TL.PPL.OlSCi.lfXLUJ.YUGfHUO 
1  ,XA,YA,GFArSP,XL 1\,YLIM,PHL IN 
CriPMON  T  (  SO)  ,\?),  SIJURCr  (  IP  )  ,  P  I  ,  K  ADf.  Cm GR  AD ,  M  1  ,NOS,  IX, 
l,XAISOl),YA(S01  l.DFArSP,  XL^NtYL  IN.r’HLIN, 


JISC.'SYDX,/, 

DC , A,YC( 501  ) 


1  PUS  r  {  SOI  )  ,  AAT  (  S('l  )  ,  TAiNfPtNsLOK  I  jNf’.LOK?  ,  Y0(  SOI  1  ,»'CX 

/k!-',  1/  ‘::>L(  )  ,XGH(  12  1,1  1. ,  S L  ,  S  ?L  ,  C 1. 2  ,  SL 2  ,  XOD  ,  YOU ,  Ri.'U  ,  C 1  , 


1  A  A  (  1  ?  )  ,  f.U>,  (  1  2  )  ,  CC  (  1  ?  )  ,  !'L)  (  1 2  )  ; !.  f-  {  1 2  )  ,  r  r  (  1  2  ) 

1,ANIA,  Ai- (SOI  )  ,cr  (  son  ,  PfM2S)  ,L  I  ,PL  (  2S)  r  NkPIPSI, 

iNr,.N'.fiS,LLy,:-'G2 


RPhLrNi  =  P'lLl\<f'i  grad 


I  -"IX 

AA{ I ) =0. 

Rh( I ) =0. 

CC( I )=0. 

DD{ 1 ) =-TAN(RPHLIN) 

FF{ I )=1. 

FF  (  I  )^-XLIN<^DD{  1  )-YLlN 

RI-.TURN 

END 


Fig.  59.  Cone  geometry  subprogram. 
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('.91'^'; 


UiUJt'L'.- 


A,  f)0,f  ,  I  (■ ,  xr.i’,  Y:/-  ,  9  i  ,  f'  ":’  ,rcc, ,  91  AC.,  f)Y:)X^ 


1  Y,  l.Sf  11  tLL  ,(.LR,  SI.  ,  ;.L;‘.C.‘i.,S2L,  1  L  t i’L  -  D  I  S',  li ,  XO' : »  Y.  IC  ,  RC:(; 

J  .XAtYA.fU  A1  S>’,;<L  r.‘,YL  ri.PHLlN 
COffn)\  T  l  ‘SOI  ,  1  in  ,  S'>U'!Cf.  (  1  ?  )  ,  PI  ,P.An;.C,nf  op  A0,P  1  ,NJf, ,  I  X,  PC  ,  A,  YC(  so  1  ) 
1  ,XA(  SCI  )  ,  YA(  ‘jOi  ),;)!■/'  1  SP,  XL  in,  YL  in,  PHL  I  M, 
i  PMS  r  (  SOI  )  ,  AAT  (  50  I  )  ,  TAPFR,NP;  DKI  ,  ,  Y0{  SO  1)  ,  ^CX 

COMMO  1  /HP,  1/  PPL  ( ;>‘n  , xr,ri(  12  j  ,tl ,Sl,spi  ,ci.?, slp, xnn,  y  in, rci'.cl, 
l  A  A  {  l  2  )  ,  on  (  1 2  )  ,  CC  I  1  ?  I  ,  no  (  1  2  )  ,  1  •:  (  1 2  )  ,  f  F  (1  2  ) 

I ,  ANf  A,  AF  «  f.oi  I  ,CF(  SOI  )  ,*-M{?S)  ,LI,PL  1  25  j  ,  Nkir  (  2*0  , 

llMR,NSiaS,LLX,Mr.? 

DIVLNSION  TI  (son  ,LG(SQH  ,YGP(  I?)  ,r>Ll  (25) 

DATA  WKI  rCR /AHksM  r/ 

DATA  SXvSf;/AIISXGP/ 
data  FNO(An/AM  NO  / 

Data  PMiN/. 0000001/ 

1  F0RnAT(9FC.6) 

2  rORPAK  16  1'^) 

4  Fnr»MAr  (  I'll  ) 

5  FOP.MAT(?X  , MINUS  IS,8X,3HWC  I5,«X,3HLL  15) 

6  FORMAT!  AX,  IIIP.OX,  1  tlX ,  OX ,  1  hY ,  dX  ,  2HT  I  ,  6X,  4i' '.FPR  ) 

7  F0RMAn2X,  ItoilX-COOROlNAJC  GEOKCTKY  POUMDAklFS  XGP  //(FIO.4)) 
ft  FORflAT!  IS,  If  10.0,  IS) 

9  FOPPA"  (  ?X  ,  2i‘.'.  F  1  0.4// (  2X,6HPL  I  FI  0.4)) 

1 1  FOfn-AriPX, /n\<e  i  is//) 

12  FORMAT  (  li-C)) 

13  FOPf'AT  (6X,  laninOK  ANGl  F  ,F10.?//) 

50  F0RMAT(4X,2H\\,HX,2Hr,r  ,BX,2HCC,yX,?H0[),8X,2HC(  ,8X,2HFF//) 

51  FORMAT (6F1 0.4//) 

52  FORMAT (2X, 7hYGP  L  FI  0.4//) 

IF{S0URCF{3).'je  .f  N'ORAO)  ••'.r,2  =  RG2  +  1 
PLI  (Mr,2)=PL{LLX) 

IF IMG2.GT.2. AOe.ASS (RLl {KG2)-PL I (XG2-1 ) ) .LT. .01 )  GO  TO  48 
fFILL X.GI  . 1 )  GO  TO  i81 
RFAL*8  RMTUI 
RMIN'J^I.D-  10 

IF(SGURCe{  3).Fr;.rN'aRAD)  GO  TO  181 
WRITfc(6,3)  NOS,MC,LT 
WRITE(6,1 1 )  (\RF{  1  )  , I^l ,MCX) 

WRITE(6, I? ) 

WRITE(6,S0) 

WRITE(6,S1  )  ( AA( I ) ,R8{ 1) , CC ( I ) , 00 ( I ) , F E  (  f  ) , F F ( I ) , I  =  1  ,  ,05) 
WRITE(6,9)A, ;PL{ I ), I=1,LT ) 

N0SS=N0S/2+I 

IFlSOURCi-' ID.F-.  ,5XG8)  GO  TO  200 
XGRd)-'^  xco 

XGB{  2  )  =Xlif:  +  nsCf;T(  Rfjn»K00-Y0(2i:Y0fJ) 

XG8{3)-  xno  , 
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DAie  -  69199 


200  CONTIfJUe 

WR  I  Ttf  (  6 , 7  )  {  XGR  (  1)  ,  I  -  I  ,  NWOS  ) 

DO  180  L=l,NOSS 
IFtBBlD.ifO.O.  )  GO  TO  182 
TLIM=DA3S{FF(L)  ) 

IF(DAHS(AA(L) l.GT.TLIMj  TLIM=  AAIL) 

IF{DABS(BB(L) l.GT.TLIK)  TLIM=  6B{L) 
ir(l)ABS(CC  (L)  )  .Gr.TLIM)  TLIM=  CC  (  L  ) 

IF(DABb(OD(L)  I.GT.TLIM)  TLIM=  ODd) 
iriDABSlEEd)  )  .GT.TLIM)  rLIM=  EEIL) 

TLIM  =  T(.lM*.OOnooi 

DISCU=( iCC(L)*XGO(L)+EE(L) )/(2.«BB(L) ) ) »« 2- ( AA { L ) « XGB { L ) »« 2+DD ( L ) 
1*XGB(L)+FF{L) 

IFIOABStDlSCLD.LT.lLIM)  DISCU=0. 

IFd.EQ.NOSS )  GO  70  68 

67  YGPtL)=-(CCd)*XGBd)+EE{Ln/(2.*BBd)  ) +OSQRT  (  DI  SCU  i 
GO  TO  183 

68  YGBd)=-{CCd)*XGBd)+EE(L)  )/(2.*BB(L)  ) -05QRT  (  D  I  JCU  i 
GO  TO  163 

182  YGfid  )=-(  AAd)»XGRd)*XGB(l.)+DDd)»XGed)  5-FFd)  )  /  ( CC  (  L  )  *XGC  (L  ) +E  E  ( 
ID  ) 

183  IF(DAfiS(YGBd)  ).l.T.l.E-3)  YGBfD  =  C. 

IFd.EC.NOSS)  GO  TO  IRO 
N1=N\G5+1-L 

YGBIM  >=-YGBd) 

180  CONTINUE 

WRITEI6,52)  (YGBI I ) ♦ I=1,NN0S) 

181  WR1TE(6,4) 

WRITEi6,13)  PLdLX) 

WRITE(6,6) 

DC  20  I -I, NR 
T2=I 

20  Y0( I )=A/(2.*ANEA)*(2.*T2-2.-ANEA) 

IFJSOURCE I  3) .FO.FNORmD)  GO  TO  48 
00  210  I=i,NR 
XA{I)=  CL»DFATSP-YO(  n*SL 
210  YA(I)=  SL*DFATSP+YU( I ) *CL 
DO  190 

DO  170  L=l,NOS 
LG(M)=L 

lF(BBd).CO.O.)  GO  TO  30 
XAC=XA(M) 

YAC=YA(M> 

AAA  =  AAd»+BBd)*TL*TL  +  CClLl»TL 

BBB=  -2.»BBd)*TL*d»XAC  +  2.»BBd  )  »rL«YAC-CCf  L)«TLeXAC  +  CC  d  )*YAC+n:) 
i  d)+E£d)*TL 

CCC--=  -2.«Pf'{  I  )*TLeYA.C*XAC-»-BBd  )-TL«TL«XACeXAC  +  BRd  )*YAC*YAC-Cfcd  )<^ 
1  TL^XAC-^EEd)  :  vac+FFIL) 
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Bbr<=RnB/(2.^fAAA) 

CCC=CCC/AAA 
DISC=RBn»BnB-CCn 
IFCDISC.LT.O.)  GO  TO  170 
X=-BBB  +  LISCRT?D1SC  ) 

IF  lX.GT.XGtiC'JOSS)  )  X  =  XGB (  NOSS  1 -RM I MN 
GO  TO  18 

30  :>-(BF(LMYO(H)/CL  +  FF(L)  )/(00(L)  +  EF{L>*TL) 

18  Y=X»TL+YC(M)/CL 

IF{XGB(L)  .LE.X,A\'n.X.Le.XG3(L+l  ).AND.YGR{L>.Gt.Y,A,'gC.Y.G£.YGB(L  +  l  ) 
DGO  TO  IB 

iF{XGB(L).GE.X.ANO.X.GE.XGB(L+n.AMD.YGB(L).GE.Y.A'JD.Y.GE.YGB(L+l  ) 
DGO  TO  15 
170  CONTINUE 

15  IF«DABS(2.«BB(L)»Y+CC(L)*X+Ee(L)  )  .LT.RMIN)  GO  TO  17 
DYDX=-(2.«AA(L)*X+CC{L)«Y  +  DD(L)  )  /  t  ?  .  *  BB  (  L  )  *  Y  «-CC  (  L  )  *X +E  E  (  L  )  ) 

TS=DYDX 

IFIOABSI  1.  >TL*TS‘i  .LT.RMIN'  GO  TO  22 
IFIY.GE.O. )GU  TO  115 

T! (M)=PI/2,-0ATAN2( (TS-Tl) » ( l.  +  TS*TL)  ) 

GO  TO  16 

115  TI (M)=PI/2.-0ATAN?( (TL-TS) , ( l.+TS*TL) ) 

16  TI (M1=DABS(RADEG*TI (M) ) 

GO  fO  21 

17  TI (M)=PL(LLX) 

GO  TO  21 

22  TnM)=0. 

21  IF(50URCE{  12). EQ, WRITER)  WRITEIG.B)  M,X,Y,TnM),L 
190  CONTINUF 

00  70  1=1. Ml 
00  70  K=1,N0S 

70  T( I ♦K)=-l. 

DO  A9  KG=liMl 

TI  {KG)  =  (  T  I  {KG)+TI  (KG+n  )/2. 

lF(lGIKG).NE.LG(K0+l),AN0.LGIKC).GT.(NnS/2))  GO  TO  71 
LG(KG)  =  ILG(KG)+l.G(KG4-n  )/2 

71  K=LG(KG) 

I  =  KG 

4E  Tn,K)=TI(KG) 

A8  RETURN 
END 
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SURROUTINC  NKJRy 

DOUR  LI:  PRECISION  A  A  ,  BR  ,  CC  »  OD  ♦  b  E  t  f-F  ,  XOB  t  YGR  ,  P  I  ,  OBri »  CCC  f  D I  SC  ,  DYDX ,  X, 
IY,TS,TI  fCL,CL2,SL,SL2,C2L,S2LtrL,KPL,DISCU,XOO,YCJO,ROO 
1  tXA.YAjDFATSPiXLIMtYLIN.PHLINtOOBST 
DOUBLE  PRECISION  TR 

COKMCN  T(  SOI ,  12  )  ,  SOURCE  (  12  )  ,  PI  »  R ADEG f  DEG'»  AD,  M I ,  NOS  f  I  X ,  MC  ,  A  ,  YC  (  50  1) 
I ,XA{ 501 ) , YA( 501 ) .DFATSP,  XL  I N , YL I N , PHL I N , DOBST , 

I  PHST(50l ) ,AATI50l) , TAPER, NBL0Kl,MRLnK2  ,Y0(50l),MCX 
COMMON  /JRM/  N(  12)  ,D(  12, 12)  ,[)CE{  12, 12)  ,T0(  12,  12)  ,FIPD{501  ),TC(501  ) 
l,FPvEQ,P0LI7 

DIMENSION  R( 12) ,01 12) ,SR( 12) ,NN( 12) ,DnD( 12) 

1  FORMAT ( F 10. 2, OX, F 1 1. 6, IX, F 10. 6, 2110) 

2  FORMAT  (bF15.6) 

3  FORMAT  ( I  15,4F15.9) 

A  FORMAT  {F15. 7,AF10.6,F15.7) 

7  FORMAT ( IHl  ) 

8  F0RM.AT(7X,  IHT,  10X,AHFIPD,9X,11HTRANS  COE  FF  ,  9X ,  I HK ,  9X  ,  I HL  ) 

999  FORMAT!'  WALL  THICKNESS  =  *F10.6/) 

DATA  PARALL/AH  PAR  / 

DATA  RMIN/. 0000001/ 

DATA  WRITER/AHWRIT/ 

DATA  FNORAO/AH  NO  / 

IF(S0URCE(3).E0.FN0RAD)  GO  TO  557 
DO  6  L=l,NnS 
NN(L)=N{L)+1 
NNL=NN( L ) 

6  nCF(NNL,L)=l. 

WRIT6(6,8) 

DO  5  L=1,N0S 
NL=N(L) 

0DD(L)=0. 
on  5  I=1,NL 

IF(L.EQ.l)  WRITF(6,999)  D(I,L) 

.  5  DnD(L)=DDD(L)+DI I ,L) 

DO  60  K^l,M,l 
DO  60  L=l,NOS 
IF  (T{K,I.)  .LT.O.  )  GO  TO  60 
IF(L.EQ.NBL0Kl.0R.L.EQ.NBL0K2)  GO  TO  100 
GO  TO  101 

100  FIPD(K)=0. 

TC(K)=0. 

GO  TO  60 

101  TH=T(K,L)*DEGRAD 
DD=2.»C0S( TH)«nnO(L)»PI 
S  =  SIN(  TH)  *SI(-M(  TH) 

SRIl  )  =  SORT(nCf:{l,L)-S) 

,  IFIPOLIZ  .EO.PARALDGO  TO  210 
RR=l SRI  1 )-COS( TH)  )/ ( SRIl )+COS( TH) ) 


Fig.  6l.  Multilayer  transmission  subprogram  -  Page  1. 
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OATr 


(.vi'Vv 


r.  n  T  n  t  i  « 

<.  A.  «. 

?1()  nR=(  SR  (  I  )-i)CK  (  l,L  I  i^COSl  Th)  )  /  (  biU  1 ) +DCl  (  t  ,L  )*Cf)S{  TH)  ) 

211  CONTINUE 
NL=N(L> 
on  10  I  =  l,iNL 
n  =  i  +  i 

S«{ I  I )  =  SORT(DCE ( I  I ,L)-S) 

G(I )=2.»PI*D( I ,L)»SR(  I  ) 

IFCPOLIZ  .EQ.PARALL )G0  TO  llO 
R( I )  =  (SRt I  I )-SRf I ) )/(SR( 11  )+SRl  I  )  ) 

GO  TO  10 

HO  R(  n  =  (OCE(  I,L)*SR(  I  I  )--OCE(  I  I  ,  L)  *SR  (  1)  )  /  (  DCE  (  I  ,  L  )  ♦S'?  (  I  I  ) +DCt  {  1  I  ,  L  ) 
1*SR(  i  )  ) 

10  CONTINUE 
AQ=1.-RR 
DO  15  I=l,NL 
15  AQ=AQ#( l.-R( I) ) 

AQ=;1./A0 

W=l. 1603i4E*4/FRt0 

GG^GIU/W 

CG=COS(GG) 

SG=SIN{G&) 

AD=PUDCEn,L)«TD(  I ,  L  )  *0  (  i  ,  L  )  /  (  W*SR  (  1)  ) 

X1=C.G*  (  1  .-  AD) 

Yl=-sr,#(  1  .-AO) 

X?=-RR#C'.t  (  1  .♦AO  J 
Y2  =  -RR*5w#  (  l.  +  AD) 

X3=-RR#CG#( l.-AO) 

Y3=RR»SG#( 1 .-AD) 

X4=CG»( l.+AO) 

Y4=SG^( l.+AD) 

NNL  =  iNN(L) 

DO  35  1=2* NNL 
IF(I-NNL)  25,20,50 
20  Ul=l. 

U2=~R(NL) 

U3=-R{NL) 

U4=l, 

V1=0. 

V2=0. 

V3=0. 

V4=0. 

GO  TO  30 
25  11=1-1 

AD=P»*[:CE{  I  ,L)»TO(  I  ,L)#D{  I,L)/(W^SR(  I  )  ) 

GG=Gl I )/W 
CG=COS(GG) 

SG=SIN(GG) 


Fig.  61.  Multilayer  transmission  subprogram  -  Page  2. 
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MKJRM 


DATE  =  6919P 


U1=CG* ( 1 .-AD J 
V1=^-SG*«  1  .-AD) 

U2=-R( I  I ) •CG*( i.  +  AO) 

V2=-R(II)*SG*^(1.+AD) 

U3=-R( I  I )«CG*( l.-AD) 

V3=R(II)*SG*(1.--AD) 

U4=CG»( l.+AD) 

V4=SG« ( l.+AD) 

30  Pl  =  Xl*Ul-Yl*\'l+X2»U3-Y2*V3 
Ql==Yl»iJl+Xl«Vl  +Y2*U3  +  X2»V3 
P2=Xl»U?-Yl*V?  +  X2*U4-Y?*V<i 
02=YUU2  +  X1*V2  +  X2*V4+Y2*U4 
P3==y3*Ul-Y3»Vl  +  X4*U3-Y4«V3 
03  =  Y3*Ul^X3<fVl  +  X4*V3+Y4*U3 
P4=X3*U2-Y3»V2+X4*U4-Y4*V4 
04=Y3*U2+X3*V2+X4*V4+Y4*U4 
Xl  =  Pl 
X2=P2 
X3=P3 
X4  =  P4 
Yl^Ol 
Y2=02 
Y3=Q3 
35  Y4=Q4 

RCR=(-X3*X4-Y3*Y4)/(X4«X4+Y4*Y4) 

RCI-(-Y3*X4+X3*Y4)/(X4*X4fY4*Y4) 

RC2  =  RCR»RCR+RCI*ftC  I 
RC  =  SQRT(f<C2) 

*»R=(Xl  +  X2ttRCR-Y2*RCI  )*AQ 
TI  =  {  Yl+Y2*RCR+X2»RCn«A0 
TC2=TR*TR+TI*TI 
TC(K )=SQRT(TC2> 

IF{TR.EQ.O,.AN-D.T' .EO.O.  )  TI  =  TI4RMIN 
XX  =  DATANI2(  TI  ,TR) 

48  FIPD(K)=-RADeG»(XX+DD/W) 

IF(FIPO(K)  .LT.O.  )  FIPD(i<)=FIPD(  K)+360. 

IF(SOURCE( 12). EQ. WRITER)  WRITE{6,1)  T { K , L  )  t F I  PD ( K ) , T G { K ) , K , L 
60  CONTINUE 
50  CONTINUE 
557  RETURN 
END 


Fig.  61.  Multilayer  transmission  subprogram  -  Page  3 
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A  T  a:>j.  u 


n  A  T  r 

i,' I  » 


/  ri  »  ^  ^ 

11  -J  I  -J  -1 


SUPRGUT  I  fd-  AIAI'PR 

D G II R.  I.  t-  P  R  h  M  S  I  n *.  A  ,  ('  ( W  C  C .  0 D ,  r  f. ,  I-  F  ,  X  G ' ' ,  Y  P. ,  i-- 1  .  f'  R I ,  C  C  C  , I  S  ,  I  ’  Y 0  A  ,  , 
I  Y,TS,  n  ,Cl.  ,CL?,SL  ,3L?,C?L,S?L,  Tl.  .KPLvDI.SCGtXUr,  YL)P,ROt; 

I  ,XA,YA,t)r  Al  bP,XLIN,YLIM,PHLlN 
CGMMO'vj  /TGR/  Y,srAP;?5) 

CO.Y.YOM  T(  ‘jOI  ,  12  )  ,  SOURCF  {  I2  )  ,Pl  ,  ff  ADF  G ,  DEGP  AO,  M I  ,\'ns,  I  X,  KC,  ^  ,  YC(  301  ) 
1  ,XA(  501  )  ,  YA  (  ‘lOn  ,  (U-ATSP,  XL  I  N ,  YL  I  N  ,  PHI.  I  , 

I  PHSri  501  )  ,  AAT(  SO  I  )  ,  T  APER ,  .NRLOK  1  ,:jyLl)K2  ,  Y0(  50  I  )  .  MCX 
11  FORMAT  {bF-lO.Sl 

152  FORMAT  (  3X,  SSHUNIFORM  APERTURE  D  I  S  TR  I  BUT  I  0\'  USED.//) 

153  FORMAT (2X,?riO. A) 

154  FORMAT (6X, ) OHORSTACLC  HEIGHT  F10.6///) 

155  FORMAT  {6X,27I!0PSTACL£  HEIGHT  MODIFIED  riO.6///) 

DATA  FLIN/AHFLIM/ 

DATA  CUS6/4HCnS6/ 

DATA  CIRCl F/6HCIRCLE/ 

DATA  CnSl/GHCOSl  / 

DATA  RAD4/6HRA04  / 

DATA  RAD5/Gh;UD5  / 

DATA  PL1/6HPL1  / 

DATA  C0S2/6MCnS?  / 

DATA  FMOAAT/OHFNQAAT/ 

DATA  FNPHST/6HFNPHST/ 

DATA  BLOCK/AH  V-LO  / 

onsT=o, 

RAD:--.5«A 

RADSQ  =  RAn#'^An 

RADI =1. /RAD 

PADI3Q=:RADI*i<ADI 

IF{SOURCE(y). HE. BLOCK)  GO  TO  13' 

ReAD(5,ll)  HOST 

WP.ITF(6, 154)  OUST 

IFISOHRCEIO)  .?1E. CIRCLE)  GO  TO  13 

RAOY=RAO/V 

OGSTSQ=ODST*OBST 

RLA0(5,U)  OBSTD 

THFT^  Y/(2.*A) 

RATI0  =  0BST*0B5.T/{RADY  +  nBSTD»SIN(  THFT  )  )»*2 
A0=  KATIO^»PI«RADSO 
THFT1=1 (V, 

THFT1=  THETUDEGRAD 
DO  15  N=1 , 10 

THET?  =  THF ( 1- ( THET 1-S I  MI THFTl )-PI * { l.-RAT 10) )/( l.-COSITHFTl ) ) 
IF{ABS(THFT?-rHFTl).LT.l.e-4)  GO  TO  17 
THFTin=  RA0EG*THL-T1 
WRITE (6, 10)  IHETID 
19  FORMAT!  'THFI A  =  'FU.O) 

15  THETi=THET? 


Fig.  62.  Aperture  taper  subprogram  -  Page  1. 
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DArC- 


69199 


17  ORST  =RAD*r.OS(  rHlfTl/;i.  '! 
wRrTf-:(6,  ib‘>)  OMsr 

ASTR  IP=(  PI-THET1+  S  IN  <  THItT  1 )  l#RADSO 
WRITE(6,18)  AQtASTRIP 

18  FORMAT  ( 'OriSTACLE  AREA  -  '  F  12 .6  f  /  •  STR  I  P  AREA  =  ’FlZ.b) 

13  no  8?  IL=lfMl 

YC( lL)=(Yn( IL)+YO( IL4 1) 1 /?. 

YCT=YC(  ID 
AnSYCT=A8S (YCT) 

IF{SOUI^CE(5).EO.FNOAAT)  CO  TO  83 
CAM=l. 

IF( SOURCE { 0) .CO.  CIRCLE )  C AM=SORT t RADSO-YCT#YCT ) 

IFtTAPER.EC.COSl  )r,0  TO  1 
IF(TAPER.E0.CnS2)r,n  TO  2 
IF(TAPFR.FQ.PLl)  CO  TO  3 
IF(TAPER.F0.RAD4)r,0  TO  A 
IFtTAPFR.EC.RAD-^ICO  TO  5 
IF(TAPER.E0.C0S6)G0  TO  6 
IF(TAPCR.EO.FLIN)GO  TO  / 

1  AAT(  IL)=S0RT(.25«A*A-YC(  ID*YC(  ID  ) 

GO  TO  84 

2  AAT(  ID=DCOS{PI»YCl  ID/A  I 
GO  TO  84 

3  IF(ARS(YC(  ID  1 .11.0.986603)  AATnD  =  l.»CAM 

IFIARSIYCI  ID  )  .G£.  0.9866  03.ANO.AbS(YC(  ID  )  .  LT  .  1 .  1  83924  )  AAT  (  ID  =  (  .9 
17B9-(ABS(YC I  ID ) -0. 98660 3 ) *0 .042  I /0 . 197 32 1 ) •CAM 
IF  (  ABSIYCI  ID  )  .GE.  1 .  I B3924  .  AMD.  AHS  I  YC  (  I  D  5  .  LT  .  1 , 38 1244  )  A  AT  I  I.D  =  (  .  9 
l368-( ABS( YC(  ID )-I.  183924 ) *0, 063 1 /O .  1 97 3 2  1  ) «C AM 
IF{  ABS(YC(  ID  )  .GE.  I  .  38  1 244 .  AMD.  AHS  {  YC  (  I  DI.LT.l. 678565  )AAT(ID  =  (.o 
17  37-{ ABSIYCI ID ) -  1 . 38  I  244 ) . 3369/0 . 1 97 32  I ) *C AM 
IFIABSIYCIID  ).GE.  1.5  7  85o5.  AMD.  AOS  (YC(  ID  ) .  LT .  1 . 7758  85  )  A  A  T  I  ID  =  (  .5 
136B-(ABS(YC( ID )-l.57H565) *0.3684/0. 197321 )*C AM 
IF( AOSIYCI ID ) . GE. 1. 775885. AMO. ACS {YC( ID ). L T . I .973206 ) AA T { IL)=( .1 
l604-( ADS ( YCl  ID )- I . 775885 ) #0.03 16/0 , 197 32  1  ) tCAM 
IF  (  ABS(YC(  ID ) .GE.  1 . 97 3206. AMD. ABS ( YC ( ID  )  .  L T . 2 . 1 70526 ) A AT ( IL)  =  ( . 1 
1368-(  ABS(  YC(  ID  )  -  I  .  97  3 .206  )  »0 .0?  lO/O  .  19  73  2  I  )*CAM 
IFIABSIYCI  ID  ).GE.2.170526.AMD.ArjS(YC(  I  D  )  .  L  T  .  2 . 76  24  68  )  A  AT  (  I  L)  =  (  .  1 
115P-{ABS(YC( ID  1-2.  170526) *0.03 16/0. 591 96?) wC AM 
If (ABSIYCI ID ) .GE.2.762408 )AATI I D -0 . 0842*CAM 
GO  TO  84 

4  AAT(  ID  =  {  l.-YCT*YCT*RADISO)*CAM 
GO  TO  84 

5  AAT(  ID  =  {  (  l.-YCT*YCT«RADISQ)**2)*CAM 
GO  TO  84 

6  AAT(  ID=DCOS(PI*YC(  I D /A  )* *2  +  5 TAP  (  MC X  ) 

GO  TO  84 

7  AAT( IL)=l-ABS( YC( ID )/RAO 
GO  TO  84 


Fig.  62.  Aperture  taper  subprogram  -  Page  2. 
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ATAPI-K 


UM  <■ 


^.,91 


«3  AAT(  ii  )  =  l. 

l(-(Aii,SYCI  .LT.lM'Sr  )  A.\r  (  11.  )^(). 

IF  t  SOURC"  (  6)  .CQ.P 'JPMST  1  GO  Tl{  B5 
PHST(  ID^kftSOP.K  1I0.?5^YC(  IL)*YC{  ID  ) 

GO  in  02 
85  PHSTI  IL  )--0. 

82  CUNTI'’JUE 

i:-(  SOURCE  (  5  )  .F0.F>J0AAT.  A^l().SOURCh  (6  )  .rO.nU’HST  )  >;P  ITL  I  6,  15  2  ) 
WRl  TF(6, 15  3)  (  AATI  IL)  ,PHS1  (  1  D  ,  I  L=  1 ,  f-U  ,  10  ) 

RETURN 

END 


Fig.  62.  Aperture  taper  subprogram  -  Page  3. 
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GRAPH 


DATE  =  69199 


SURROUTINF  GRAPH { X , XMAX , XM I N , NUMIFR , R , I  SC AL E J 
DirFNSIDH  X(.MUNBtR)  ,n(?6)  ,P{4)  ,  S  (  1  ) ,  B I  NUFlihR  ) 

DATA  G/'  +•/ 

DATA  S{  I  )  ,P(  1 )  ,P(2)  ,Pl  3)  ,P{4)/ '  *,‘*'t*  **»*  *•/ 

IF( ISCALE.EO.O)  ISCALE=l 
FACTOP.=  lOO./(  XMAX-XMiN) 

WRITE(6,31) 

31  FORMAT { IHI//20X, 'FAR  ZONE  POWER  PATTERN  (D0)< 

DO  2  K=l,?6 

2  0{K)=S(l) 

L=0 

DO  5  N=l, NUMBER 

63  L=L+1 

IF(MOO(L-l,ISCALE) )  50,50,51 
51  WRITfc(6,53) 

53  FORMAT! IH  ) 

GO  TO  55 
50  J=1 
K=1 

Y=(X(N)-XMIN)«FACTOR 
IFIY.EQ.O.)  Y=1E-10 
Y=  20.«AL0G10{ Y)+60. 

10  IF(Y.LT.3.5)  GO  TO  20 
Y=Y-4.0 
J=J+l 
GO  TO  10 

20  1F(Y.LT.0.5)  GO  TO  30 
Y=Y~l.O 
K^K+1 
GO  TO  20 
30  0(J)=P(K) 

WRITc{6,3?)  R{N),0 

32  FORMAT(3.<,F10.4,2X,26A4) 

0( J)=S( 1) 

55  IFCMOnCL-l, 10) 140,40,64 
40  WRITE (6, 42)  (G, I-l, 11) 

42  FORMAT! 1M+,6X, IIAIO) 

64  IF(M0D!L-1,ISCALE) )  5,5,63 
5  CONTINUE 

RETURN 

END 


Fig.  63.  Graph  subprogram . 
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SLi 


0  ^  I  I 


J') 


'ill  Si:  (Si  ,  XMA/.vijr-r.i-'^ ,  ■, ,lSLMS>’,SL^',n, 1  ) 

PI  f-T  \  S  K'  sr  (  •  U  ,  A  !  'HJ-hf  P  )  ,  PLi'  (  r )  ,  I'  (  PU'  W  ''  ) 

P  1  =  3.  I'.  1  \2 

01  GR  A(.i^P  I  /  l(;U. 

RAI,'1-G^U10./P1 
11  Ff'PMAT 

PAT  A  ^'^•^7(J/ 

;  (RCf-'.co.  1 )  GG  rn  12 
0 

K  =  0 

I=KU‘-'r,Fr’:-l 
DC)  1  :nI=2,I 

IP(Ac(G-i  ),  LI.  AL  A\iD.AF(Mn).Ll  .AF(\')  )  '7)  TO  2 
GO  TO  1 
?  M  =  M+1 

A{M)'-(AL(\')/Xf-'AX)*^100. 

SLG( V)=R(M) 

writ: (6, 1 1  )  AC  M , SLnCK) 

1  CCiMTI.GUt 

ir-CMOl.FO.  l.AG'n.R.GE.Gi  GO  TO  >3 
GO  TO  14 

13 

Sim  --SLP  ( f'l-  ) 

WRIT^-Co,.?'))  •••'01 

25  rOU-ATC  .”r;i  =  *15,'  M.V- .M='15//) 

14  sl:'ax.^o.o 
DO  3  M=l,'l 

1F(A(;-!).GT.SL‘'1AX.  AOD.ACNl.LT.T^.  )  GO  TO  G 
GO  TO  3 

8  SLOAX=A(N) 

PSLf7,X-SLC'.  (0)  *OfcGRAD 
WRITK65II)  r.SLOAX 

3  CONTI GUF 

irCSLMAX.LT.O.OOOnOOn  go  TC.  6 
5108  =  20.  V  ^  logic  (  I'SO./SLOAX) 

GO  TO  7 
6  WR  IT  F  { 6 , 9  ) 

9  FORMAT!  •  S'O  GIDFLOGF  DETCCTED' 

HSLl=5Lri«''c'.RAD 
RSLNAX-=HSL  1 

Rt-TURO  2 

7  WRPL(6,5)  SL  MAX,  SLOB 

5  FORMAT  (  lOX  ,  2vH  f'AXlMUf-'  SIOCLOP.L  LEVEL  ,  F  i  0. 4 , 5  X  ,  BMP  L  RCEO  f  .  ,  vX  ,  I H 


1  ,F10 

.4, 

6H  Dh  ) 

P,VW  = 

1 

RE  TO 

12 

Vr/f/r 

0 

P51. 1 

--■SLBiniEO'-AO 

BMP 

BSL  1 

50 

Rf  lU 

’n 

1 

END 

Fig.  64.  Sidelobe  subprogram .. 
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HPl^W 


DATE  =  0919'^) 


SUPRDUT  IME  HPr>  W(  AC  ,  y>'AX,''lUMfiCri ,  B) 

D I  1  CE>)  AE  ( 'lUyBC  P  )  ,  C-  (  >JUNBCR  ) 

?F  (  AFJ  N)  !ge!'!k  707  UXMAX.  AND.  AE  (  N-  1)  .LT  .  0  .  ^  xUJ!!  J 

IF(  AE(N)  .LT.0.707UXMAX.  ANn.AE(N-l  )  .oE.0.70.  l«XMAX  J 

GO  TO  1 

Al=Ab{N) 

A2=AriN-l ) 

PA1=RIN) 

PA2=B(N-1 ) 


GO  TO 
GO  rn 


10 

20 


GO  TO  2 
A3=AE(N) 
A4=At ( N-1 ) 
PA3=fi(N) 
PAA=B( N-l » 


COf'lTU^UE 

AN'.Gl  =  PA2-  {  PA2-PAI  ) 
At'3G2  =  PA3+(?A4-PA3) 

bw-angi-ang? 

WR1TF16,40)  b’.V 
FORE'.AT  (  ///  lOX,  25H 


»l,7071*XMAX-A?)/t  M-A25 
*(  .?071*XKAX-A3)/{A4-A3) 

HALF' POWER  t^FAMWlDTH  ,F10,4, 


3X, 


BHOEGREFS.//) 


RETUK^J 

END 


Fig-  Kali-.'-'-wer  bc&mwidth  subprogram. 


I 

I 


X  X 


c 

'3- 


•> 

AM 

■Mk 

c 

•> 

AM. 

m 

.k» 

in. 

II 

o 

*m| 

c 

Lm 

o 

a 

>c 

o 

■«^ 

UJ 

o 

5: 

— >• 

Ii 

MM 

to 

M 

mJ 

H 

o 

to 

J, 

•» 

lO 

LL 

Uj 

UJ 

• 

to 

KJ 

w* 

rM« 

*u 

X 

o 

XJ 

uJ 

< 

l-M 

:i 

MM 

■X 

MM 

MM 

s: 

a 

to 

c« 

*. 

M 

lO 

tO 

«> 

X 

u* 

.J 

y 

X 

y 

t- 

Q 

— .k 

-J 

>• 

‘JU 

.u 

LO 

•> 

C 

o 

X 

O  2^ 

a: 

w 

o: 

•i 

X 

o 

UJ 

•> 

• 

w 

AM 

A 

il 

—A 

M 

a. 

<I 

LU 

II 

l< 

CL 

UJ 

<u 

^4 

II 

II 

u 

3. 

U. 

a 

nj 

o: 

to 

3j 

•» 

CL 

UJ 

ft 

, 

AM 

1-5 

to 

AM 

MM 

.J 

X 

CL 

c-O 

^s. 

« 

II 

h- 

<•* 

MM 

l; 

C’ 

Cl 

|i 

<■ 

J. 

•> 

r' 

o 

_J 

M 

X 

AA. 

c. 

iTi 

U- 

* 

to 

M 

L> 

;> 

AM 

X 

X. 

•> 

o 

LL 

UJ 

AM 

M-* 

o 

a4 

#«* 

••4 

UJ 

*X’ 

nr 

f\ 

C.' 

AM 

mJ 

M 

li¬ 

S' 

o 

a: 

w 

O’ 

or 

X 

C 

n* 

—* 

-M- 

» 

AM 

ar 

m 

a: 

#> 

#• 

*■ 

to 

tO 

c 

•k 

MM 

> 

AA 

M 

w 

*>■ 

c 

.J 

o 

« 

M 

>- 

r'^ 

<7 

¥ 

t/. 

rW 

M 

II 

<o 

V 

f  J 

ii 

a. 

VO 

II 

/■^ 

nj 

• 

-mJ 

fvi 

mJ 

fL 

M 

.M 

to 

It 

mJ 

u. 

w 

a 

> 

4-4 

> 

•» 

II 

•i. 

•M« 

h- 

X 

II 

♦ 

s 

o 

II 

o 

«• 

O 

«J 

o 

II 

•J 

y* 

Cj 

MM 

MM 

mJ 

u 

Ma 

m’ 

AM 

• 

>• 

-J 

t. 

V' 

*4 

t  • 

> 

• 

It 

II 

s 

u 

or 

tJ 

LL' 

9 

o 

cr 

zz 

y 

O 

o 

MM 

u 

J.) 

lU 

a' 

•iH* 

-J 

<r 

X 

Ji 

uo 

MM 

M 

aJ 

X* 

AM 

a 

o 

ar 

L  J 

< 

M 

•> 

A 

-J 

a>4 

•M 

a 

q: 

il 

IL' 

<r 

«j 

< 

CL 

CD 

to 

< 

M 

M 

MA 

X 

1— 

o 

a 

AM 

a 

♦ 

•k 

u 

U. 

A 

II 

CL 

U' 

a" 

AM 

Ii 

o 

1/“' 

o 

'kO 

it 

to 

O'' 

o 

t< 

u: 

X 

•> 

C 

<3 

II 

*7^ 

<r 

♦ 

4 

M 

•<** 

•• 

p. 

a 

>z 

11 

L' 

h- 

• 

ii. 

X 

tJ 

III 

MM 

o 

La 

u 

9 

<t 

to. 

UL 

y 

X 

t 

<MM 

u. 

c. 

\/^ 

ti 

lO 

i// 

Ii 

c. 

ti 

n 

*• 

u 

»- 

cr 

c 

M-4 

o 

A 

•> 

.  J 

X 

>- 

M 

X 

•* 

f- 

< 

ui 

w 

< 

\L 

vO 

c 

LU 

< 

<2 

i;r 

t/> 

to. 

ID 

X 

•• 

tr 

ctr 

a 

of 

to 

Ii 

P 

a 

> 

CO 

O' 

II 

Q 

t 

1! 

c 

II 

c: 

»-4 

JJ 

a 

to 

w 

> 

z> 

u. 

to 

1! 

w 

Lm 

to 

o 

LL 

t/") 

X 

t/ 

P' 

u. 

u. 

> 

II 

;o 

o 

CL 

II 

a 

II 

> 

t'. 

X 

1’ 

> 

Ia 

X 

LU 

> 

It 

> 

to 

X 

c 

II 

A 

u; 

UJ 

X 

O 

V) 

o 

2iC 

4/0 

w-^ 

y. 

X 

c 

o 

LL 

ii 

o 

t^ 

o 

«— 

AM 

5l' 

<r 

o 

to 

il 

<r’ 

tO 

Cl 

o 

o 

►— 

o 

/*> 

>■ 

V 

c 

o 

M 

< 

;r' 

e 

> 

K 

c 

9 

c 

C 

X 

rv- 

M 

e" 

M 

to 

c 

in 

</0 

i/*) 

MM 

CkT 

u 

c: 

r“* 

M— 

vo 

i'D 

>» 

CD 

to 

Cv. 

2: 

a 

r- 

•• 

O  ‘/I 

it 

C 

to 

o 

X 

o 

Ct 

a 

o 

C 

b 

;r 

51 

c 

C 

Cl 

Q 

•r 

< 

u 

a 

AM 

o 

w 

a 

a 

Mm 

c 

c- 

o 

I_ 

c 

•> 

UJ 

u: 

-L 

It 

o 

c 

r  'k 
4.M 

o 

mJ 

a 

o 

o 

o 

X 

u 

cz 

a 

1- 

*— 1 

X 

C 

UJ 

o 

>- 

O  to 

mM 

X 

»a"4 

-A 

rv, 

•> 

U) 

jj 

to 

to 

LT 

X 

o 

4 

Ni. 

L5 

>f) 

c:? 

MM 

O' 

o 

X 

> 

> 

f 

u. 

CO 

C.' 

MM 

•-M 

vO 

‘X 

o 

o 

c 

uO 

lO 

o 

—) 

c 

*— 4 

to 

cc 

t 

f\J 

U' 

MM 

2. 

Of 

t 

X 

u. 

1— 

u 

X 

X 

c 

o 

Cl 

« 

Ca 

• 

n. 

.•r 

c. 

mrni 

mJ 

Ci 

LU 

u 

'J 

a') 

VO 

fO 

a 

cc 

LL' 

O 

\/i 

C. 

40 

c 

LL 

•— • 

a. 

to 

to 

to 

UJ 

o 

• 

o 

• 

« 

»- 

J: 

> 

5 

> 

u. 

X' 

X 

> 

> 

> 

y 

o 

1- 

C' 

h- 

c 

ni 

X 

i/’t 

O 

to 

O 

to 

U' 

t/. 

-J 

(O 

lO 

to 

V*) 

.J 

U- 

c 

X 

c 

o 

V 

V, 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

>V 

X 

X 

X 

X 

X 

X 

X 

v». 

■' 

'*•  - 

s,^ 

'v 

X 

X 

X 

X 

X 

> 

X 

X 

X 

X 

X 

•- 

X 

X 

C 

O 

CS 

r-t 

a 

s 

o 

u 

hK 

IS 

W 

X 

X 

D 

to 

O 

t-i 

o 


O 

•-J 

•w' 

at 

nf 

p* 

ao 

d 

.4 

t— ! 

o 

d 

o 

U 

x> 

o 

*-3 

H 

OJ 

d 


a> 

d 

T3 

U 

U 


3 

cr 

0)  h 

oi  0, 


•.D 


ao 


il2 


on  IBM  SYSTEM  360/75. 


REFERENCES 


1.  Van  Doeren,  R.E.,  "Application  of  an  Integral  Equation  Method 
to  Scattering  from  Dielectric  Rings,  "  Report  2485-2,  19  April 
1968,  ElectroScience  Laboratory,  The  Ohio  State  University 
Research  Foundation;  prepared  under  Contract  N62269  ^T-C-OjSZ 
for  the  Department  of  the  Nav>',  U.S.  Naval  Air  Development 
Center,  Washington,  D.C. 

2.  Hahn,  G.M..  and  Pinney,  E.J.,  "Three  Dimensional  Error 
Prediction,  "  Proceedings  of  che  OSU-WADC  Radome  Symposium, 

The  Ohio  State  University,  (August  1956). 

3.  Kilcoyne,  N.R.,  "An  Approximate  Calculation  of  Radome  Bore- 
sight  Error,  "  Presented  at  the  Ninth  Electromagnetic  Windov/S 
Symposium,  sponsored  by  the  U..S.  Air  Force  Avionics  Laboratory, 
and  The  Georgia  Institute  of  Technology,  Atlanta,  Georgia, 

(14  June  1968). 

4.  Richmond,  J.H.,  "Efficient  Recursive  Solutions  for  Plane  and 
Cylindrical  Multilayers,"  Repo-  t  1968-1,  10  August  1965,  Electro- 
Science  Laboratory,  The  Ohio  .‘  tate  University  Research  Foundation; 
prepared  under  Contract  NOw  <  5-0329-d  for  Bureau  of  Naval  Weapons, 
Washington,  D.C.  (AD  624  19^). 

5.  Styron,  J.B.,  and  Hoots,  L.C.,  "Antenna  Distributions  in  High 
Perfcrmance  Aircraft  Radomes,  "  presented  at  the  Fifteenth  Annual 
Symposium  on  USAF  Antenna  Research  and  Development,  University 
of  Illinois,  (October  1966). 

6.  Collier,  J.R.,  "Effects  of  Antenna  Aperture  Blocking,  " 

Report  1180-9,  30  September  1962,  ElectroScience  Laboratory, 

The  Ohio  State  University  Research  Foundation,  prepared  under 
Contract  AF  33(616)-7614  for  Aercmautical  Systems  Division, 

AFSC,  U.S.  Air  Force,  Wright-Patterson  Air  Force  Base,  Ohio. 

7.  Breeden,  K.H.,  et  al.,  "Millimeter  Radome  Design  Techniques,  " 
Report  Number  G. I.  T.  68-38,  February  1968,  Engineering 
Experiment  Station,  Georgia  Institute  of  Technology. 

8„  Corning  Glass  Works  Materials  Specification,  Section  2,  Sheet  la, 
Material  No.  9606,  (March  9,  1961). 


113 


9.  Person?.]  communication  from  W.  Beamer,  U.S.  Naval  Air 
Development  Center.  Johnsville,  Warminster,  Pennsylvania, 

(22  September  196&). 

10.  Personal  communication  from  W.  Beamer,  U.S.  Naval  Air 
Development  Center,  Jobnsville,  Warminster,  Pennsylvania, 

(19  June  1969). 

11.  "Electrical  Test  Report  for  the  U.S.  Naval  Air  Development 
C'nter  High  Temperature  Radomes.  "  Brunsv/ick  Report  BR-794- 
124-003,  9  August  1966,  Brunswick  Corporation  Defense  Products 
Division,  Marion,  Virginia. 

12.  Personal  communication  from  W.  Beamer,  U.S.  Naval  Air 
Development  Center,  Johnsville,  Warminster,  Pennsylvania, 

(20  January  1969). 


114 


UNCLASSIFIED 

Security  Classificaticn 


DOCUMENT  CONTROL  DATA  >  R&D 

fSecurtly  ciossifieation  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  is  classified) 


I.  ORIGINATING  ACTIVITY  ^Cerpofote  outAot^Jiiiectrobcience  Daboratorv  20.  report  security  classification 

Department  of  Electrical  Engineering,  The  Ohio  | _ Unclassified 

State  University,  Columbus,  Ohio 


J.  REPORT  TITLE 

A  Two-Dimensional  Ray- Tracing  Method  For  The  Calculation  of  Radome 
Boresight  Error  And  Antenna  Pattern  Distortion 


4,  DESCRIPTIVE  NOTES  (Typt  of  report  and  mclustue  dales) 

Technical  Renort 


S.  AUT  HO  Pi'S;  (Last  names  first  name,  intliol) 


Kilcoyne,  N.R. 


f.  REPORT  DATE 

2  October  1969 


so.  CONTRACT  OR  GRANT  NO. 


Contract  N00019-69-C-0325 

b.  PROJECT  NO. 


Ta  TOTAL  NO.  OF  PACES  TtIT"*  '  ~jf  REFS 

114  12 


(a  ORIGINATOR'S  REPORT  NUMBERTS; 

ElectroScience  Laboratory  2767-2 


*b  OTHERREPORT  Ha'S)  (Any  Other  numbers  that  muy  be 
assigned  this  report) 


It.  SUPPLEMENTARY  NvfcS  It  SPOl.VJRINC  MILI1  ARY  ACTIVITY 

Department  of  the  Navy 
Naval  Air  Systems  Command 

_ Washington.  D. C_.  20360 

>J.  ABSTRACT 

A  tv'O-dimensional  ray  tracing  analysis  for  the  calculation  of  radome 
boresight  error  and  antenna  pattern  distortion  is  presented  here.  Emphasis 
has  been  placed  on  the  development  of  a  method  having  considerable  flex¬ 
ibility,  so  as  to  enable  application  of  the  method  to  a  wide  range  of  antenna- 
radome  problems,  and  on  relative  ease  of  calculations,  so  as  to  minimize 
calculation  time.  Several  example  problems  are  calculated  to  demonstrate 

the  usefulness  of  the  approach.  Comparisons  between  calculations  and 

\ 

measurements  have  been  included  whenever  measured  data  .ere  available. 
Instructions  for  use  of  this  completely  computerized  method  are  included 
along  with  several  tables  describing  variables  and  the  complete  computer 
program  with  necessary  subroutines.  Programs  are  written  in  Fortran  IV 
language  suitable  for  use  on  the  OSU  version  of  the  IBM  system  360/75  (som 
minor  changes  may  be  required  for  use  on  other  360/75  installations). 


DO 

I  JAN  *« 


UNCLASSIFIED 

Security  Classification 


UNCLASSIFIED 

Security  Classification 


Two-- dimensional 
Ray-tracing 
Boresight  error 
Pattern  distortion 
Source  taper 
Aperture  blocking 
Hyper- environment 
Design 
Analysis 

Computer  program 


INSTRUCTIONS 


1.  ORIGINATING  ACTIVITY:  Enter  the  name  and  address 
oi  the  contractor,  aubcontraetcr,  grantee,  Department  of 
Defense  activity  or  other  organization  (corporaie  author) 
issuing  the  report. 

2o.  REPORT  SECURITY  CLASSIFICATION:  Enter  rhe  over¬ 
all  security  classification  of  the  report.  Indicate  whether 
"Restricted  Data"  is  included.  Marking  is  to  b.-'  in  accord¬ 
ance  with  appropriate  security  regulation:. 

26.  GROUP:  Automatic  downgrading  is  specified  in  DoD 
Directive  S200. 10  and  Armed  Forces  Industrial  .Manual. 

Enter  the  group  number.  Also,  when  applicable,  show  that 
optional  markings  have  been  used  for  Group  3  and  Croup  4 
as  authorized. 

3.  REPORT  TITLE:  Enter  the  complete  report  tttle  in  all 
capital  letters.  Titles  in  all  cases  should  be  unclassified. 

If  a  meaningful  title  cannot  be  selected  without  classifica¬ 
tion,  show  title  cli'ssification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DEISCRIPTIVE  NOTES:  If  appropriate,  enter  the  type  of 
report,  e.g.,  iiUerit.-i,  progress,  summary,  annual,  or  final. 

Civs  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AITHOR(S):  Entci  the  namefs'  of  aulh-.i-fs)  as  shown  on 
or  in  the  report.  Enter  last  name,  first  name,  middle  initial. 

If  military,  show  Mnk  and  branch  of  service.  The  name  of 
the  principal  sutho'  .s  an  absolute  minimum  requirement. 

6.  REPORT  DATE:  Enter  the  dale  of  the  report  as  day, 
month,  year,  or  oirnth,  year,  if  more  than  one  date  appears 
on  the  report,  use  date  cf  publication. 

7a.  TOTAL  NUMHER  OF  PACES:  The  total  page  count 
should  follow  normal  pagination  procedures,  i.e.,  enter  the 
number  of  pages  containing  information. 

76.  NUMBER  OF  RFFERENCES:  Enter  ihe  total  number  of 
references  cited  in  the  report. 

8o.  CONTRACT  OR  GRANT  NUMBER:  If  appropriate,  enter 
the  applicable  numbei  of  the  contract  or  grart  under  which 
the  report  was  wrhten. 

86,  8c,  £  8d,  FROJFXT  NUMBER-  Enter  the  appropriate 
military  department  i-lcntification,  such  as  project  number, 
aubprojcct  number,  system  numbers,  task  number,  etc. 

9a.  ORICIXATOR'-S  REPORT  MMBERiS)-  Enter  the  offi- 
cial  report  number  b-e  which  ihc  document  will  be  identified 
and  controlledb)  the  originating  activity.  This  number  must 
be  unique  to  this  report. 

96.  OTHER  REPORT  .NUMBF.RfSI-  If  the  report  has  oren 
assigned  any  other  -t-port  numbers  (either  by  ihr  originator 
or  by  the  sponMor),  a, so  e-.ircr  this  njmbcrfsi. 


10.  AVAILABILITY, 'LIMITATION  NOTICES:  Enter  any  limi¬ 
tations  on  fuithet  dissemination  of  the  report,  other  than  those 
imposed  by  srcurily  classilicalion,  using  standard  slatemcnta 
such  as; 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC.” 

(2)  “Foreign  announcement  and  dissemination  of  th's 
report  by  DDC  is  not  authorized,** 

f3)  *‘U.  St  Government  agencies  may  obtain  copies  of 
this  report  directly  Pom  DDC.  Other  qualified  DDC 
users  shall  request  through 

(4)  *‘U.  S.  military  agencies  may  obtarn  copies  of  this 
directly  from  0DC»  Other  qualiiicd  users 
shall  request  through 

_  _  »• 

^5)  distribution  of  this  report  is  controlled,  Quali- 
tied  DDC  users  sholl  request  through 

•t 

H  the  report  hns  been  furnished  to  the  Office  of  Technical 
Services,  pepartnvsnl  of  Commerce,  for  sale  to  the  public,  indi¬ 
cate  this  laci  anc  enter  the  price,  if  known. 

11.  SUPPL.EMF.?iTARY  NOTES:  Use  for  addi'ional  explana¬ 
tory  notes. 

12.  SPONSOF.INO  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmen  al  project  office  or  laboratory  sponsoring  (pay- 
ing  for)  the  reietich  and  development.  Include  addrc.ss, 

?3,  ABSTRACT:  Enter  in  abstract  giving  a  brief  and  factual 

su.mmary  of  vh  »  documeni  indicative  of  the  report,  even 
thoueh  it  mav  also  appear  elsewhere  in  the  body  of  the  tech¬ 
nical  fcp*'’*.  If  additional  space  is  required,  a  conlinoaiion 
sheet  shall  le  attached. 

It  is  iiighly  desirable  that  the  abstract  of  classified  re¬ 
ports  be  unclassified.  Each  paragT»iph  of  the  abstract  shall 
end  with  an  indication  of  the  military  security  classification 
of  the  information  in  ihe  paragraph,  represented  as  (TS).  (S) 
(C),  or  (U). 

There  ic  no  limitation  on  the  leneth  of  the  abstract.  How¬ 
ever.  the  suggested  length  is  from  1^  to  22.‘>  words. 

IL  KEYWORDS:  Key  words  are  technicalW  meaningful  terms 
or  short  phrases  lha:  characterize  a  repori  and  may  be  used  as 
index  entries  for  cataloging  iKc  report.  Key  \^ords  must  be 
selected  so  that  no  security  classification  is  required.  Identi¬ 
fiers,  such  as  equipment  model  dcsicnaiion,  i.*ade  name,  null- 
t^ary  project  rode  name,  ‘ogranhic  focalion,  may  be  used  as 
key  weras  hut  will  be  followed  by  an  indication  of  technical 
context.  'Ihe  assignment  of  l.nks,  rules,  and  nelghts  is 

ODtinral, 


